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Abstract. Plasma-boundaries floating in an ionized gas are usually negatively charged. They accumulate 
electrons more efficiently than ions leading to the formation of a quasi-stationary electron film at the 
boundaries. We propose to interpret the build-up of surface charges at inert plasma boundaries, where 
other surface modifications, for instance, implantation of particles and reconstruction or destruction of the 
surface due to impact of high energy particles can be neglected, as a physisorption process in front of the 
wall. The electron sticking coefficient s e and the electron desorption time r e , which play an important role 
in determining the quasi-stationary surface charge, and about which little is empirically and theoretically 
known, can then be calculated from microscopic models for the electron- wall interaction. Irrespective of 
the sophistication of the models, the static part of the electron-wall interaction determines the binding 
energy of the electron, whereas inelastic processes at the wall determine s e and r e . As an illustration, we 
calculate s e and r e for a metal, using the simplest model in which the static part of the electron-metal 
interaction is approximated by the classical image potential. Assuming electrons from the plasma to loose 
CZ3 ■ (gain) ener gy at the surface by creating (annihilating) electron-hole pairs in the metal, which is treated 

as a jellium half-space with an infinitely high workfunction, we obtain s e ~ 10~ 4 and r e « 10~ 2 s. The 
i-C , product s e T e ~ 10 -6 s has the order of magnitude expected from our earlier results for the charge of dust 

particles in a plasma but individually s e is unexpectedly small and r e is somewhat large. The former is 
a consequence of the small matrix elements occurring in the simple model while the latter is due to the 
^s^j , large binding energy of the electron. More sophisticated theoretical investigations, but also experimental 

^ ■ support, are clearly needed because if s e is indeed as small as our exploratory calculation suggests, it would 

' have severe consequences for the understanding of the formation of surface charges at plasma boundaries. 

To identify what we believe are key issues of the electronic microphysics at inert plasma boundaries and 
0^ ' to inspire other groups to join us on our journey is the purpose of this colloquial presentation. 



PACS. 52.27.Lw Dusty or complex plasmas - 52.40.Hf Plasma-material interaction, boundary layer effects 
■ - 68.43.-h Chemi-/Physisorption: adsorbates on surfaces - 73.20.-r Electron states at surfaces and interfaces 

On ' 
O ■ 

^ | 1 Introduction chemical reactor all constituents are neutral, a gas dis- 

• i-h . charge contains also charged constituents. There are thus 

^ i Low-temperature plasma physics is undoubtedly an ap- at least tw0 macroscopic scales: the electromagnetic scale, 

& Plied science driven by the ever increasing demand for where screening and sheath formation takes place [HE], 

plasma-assisted surface modification processes and envi- and the extension of the vessel. Since the observed physical 

ronmentally save, low-power consuming lighting devices. Properties of a gas discharges emerge from processes oc- 

At the same time, however, the physics of gas discharges curnn S on at least three different length (and time) scales 

is rich on fundamental problems which are of broader in- " one microscopic and two macroscopic scales - the start- 

t eres t ing point of any quantitative description is a multiple-scale 

t— i r i r j- i analysis even if it is not explicitly performed. Being exter- 

trom a lormal point of view, a gas discharge is an , . r i r 

hj. , ,j ,. ,,. , nally driven, low-temperature plasmas are moreover tar- 

externally driven bounded reactive multicomponent sys- _ , , ' . , ,., ,. . 

, T , , . i ., , , j . , • ,i oft thermal equilibrium and like other dissipative systems 

tern. It contains, besides electrons and ions, chemically ^ ... rm 

, . , , / 11+ i • j -,i feature a great variety oi sell-orgamzation phenomena kM 

reactive atoms and/or molecules strongly interacting with _ _,. „*= , ,. , r r , . • 

u ,t , J i w ,. , ,n 4 . finally, and this sets the theme of this colloquium, 

each other and with external (wall oi the discharge vessel) 7* J: ., . . , 

n . , w , • j r j j.- 1 \ u j low-temperature gas discharges, m contrast to magnet- 

as well as internal (ran to urn-sized solid particles) bound- . „ n , , • , , r <• ■ 1 i- 

T 1 ■ ,1 , n; • ically conhned high-temperature fusion plasmas, are di- 

aries. Like m any reactive system elementary collision pro- 111 • • 1 • mi 

/,,•.,,. j . \ • • rectly bounded by massive macroscopic obiects. thus, they 

cesses (elastic, inelastic, and reactive), occurring on a mi- J , . ., 

1 j , , 1 strongly interact with solids, 

croscopic scale, determine m conjunction with external J 

control parameters the global properties of the system on The plasma-solid interaction is of course at the core 

the macroscopic scale. However, whereas in an ordinary of all plasma-assisted surface processes (deposition, im- 



2 



F. X. Bronold et al.: Physisorption kinetics of electrons at plasma boundaries 



plantation, sputtering, etching, etc.) [5]. Of more funda- 
mental interest, however, is the situation of a chemically 
inert (i.e., no surface modification due to chemical pro- 
cesses, no reconstruction or destruction of the surface due 
to high-energy particles etc.) floating surface, where the 
interaction with the plasma leads only to the build-up of 
surface charges and thus to a quasi-two-dimensional elec- 
tron film which may have unique properties similar to elec- 
trons trapped on a liquid helium surface [6j or to electrons 
confined in a semiconductor heterojunction [7J. 

In plasma-physical settings surface charges play a role 
in atmospheric plasmas, where the charge of nm-sized 
aerosols [8] is of interest, in space bound plasmas, where 
surface charges of spacecrafts [9lfl0] and of interplanetary 
and interstellar dust particles [TTUT2"] have been exten- 
sively studied, and in laboratory dusty plasmas, where 
the study of self-organization of highly negatively charged, 
strongly interacting /im-sized dust particles became an ex- 
tremely active area of current plasma research [13,14,15, 
16, 17, 18, 19 . Surface charges affect also the physics of di- 
electric barrier discharges - a discharge type of huge tech- 
nological impact [SolEBEaEaEllES]. 

That surface charges at plasma boundaries could be 
considered as a thin film of adsorbed electrons ( "surface 
plasma") in contact with the bulk plasma was originally 
suggested by Emeleus and Coulter in connection with their 
investigations of wall recombination in the positive col- 
umn • Later, Behnke and coworkers [37] used this idea 
to phenomenologically construct boundary conditions for 
the kinetic equations describing glow discharges and Ker- 
sten et al. [28) employed the notion of a surface plasma to 
study the charging of dust particles in a plasma. 

Although the surface plasma as a physical entity with 
its own physical properties is implicitly contained in these 
investigations, a microscopic description of its formation, 
dynamics, and structure was not attempted. First steps 
in this direction were taken by us in a short note [29] . 
The purpose of this colloquium is, on the one hand, to 
extend these considerations, in particular, to identify the 
surface physics which needs to be resolved before a quan- 
titative microscopic theory of the surface plasma can be 
constructed and to convey, on the other hand, our convic- 
tion that the concept itself is not empty. On the contrary, 
it puts questions center stage which are of fundamental 
interest. To list just a few: 

• What forces bind electrons and ions to the plasma 
boundary? 

• How do electrons and ions dissipate energy when 
approaching the boundary? 

• What is the probability with which an electron sticks 
at or desorbs from the boundary? 

• What is the density and temperature of the surface 
plasma and are there any collective properties? 

• What is the mobility for the lateral motion of 
electrons and ions along the wall and can it be 
externally controlled? 

• How does all this affect electron-ion recombination 
and secondary electron emission on chemically inert 
plasma boundaries? 




Fig. 1. Illustration of the elementary surface processes leading 
to the build-up of a quasi-stationary surface plasma at an inert 
plasma boundary. 



The elementary processes responsible for the formation 
of a surface plasma at an inert plasma boundary are shown 
in Fig[T] Electrons and ions are collected from the plasma 
with collection fluxes jg° n = s ej ijg I f sma , where s e ^ are the 

sticking coefficients and j^f™ 3, are the fluxes of plasma 
electrons and ions hitting the boundary. Electrons and 
ions may thermally desorb from the boundary with rates 
T~l, where are the desorption times. They may also 
move along the surface with mobilities ^ e ,i, which in turn 
may affect the probability with which ions recombinc 
with electrons at the wall. All these processes occur in 
a layer whose thickness d is at most a few microns, that 
is, on a scale where the standard kinetic description of 
the gas discharge based on the Boltzmann-Poisson system 
breaks down. Thus, the above listed questions can be only 
addressed from a quantum-mechanical point of view. 

Of particular importance for the quantitative descrip- 
tion of the build-up of a surface plasma are the stick- 
ing coefficients s e ^ and the desorption times T e j. Little 
is quantitatively known about these parameters, in par- 
ticular, with respect to the electrons. Very often, s e ~ 
Si pa 0.1 — 1 and r" 1 = t+ = is used without fur- 
ther justification. Below, we sketch a quantum-kinetic ap- 
proach to calculate s e and r e from a simple microscopic 
model for the plasma boundary interaction which treats 
the interaction of electrons with plasma boundaries as a 
physisorption process [^[^[32U3MMlJ^[36l[37l[38U3M^ 
in the polarization-induced attractive part of the surface 
potential. Electron surface states ^1I4^|4^1|4^1|431H6W1 
48, 49, 50, 51, 52, 53, 54, 55 a , at most a few nm away from 
the boundary, will thus play a central role as will surface- 
bound scattering processes which control electron energy 
relaxation at the surface and thus electron sticking and 
desorption. 

Although the forces and scales are different for ions, 
they behave conceptually very similar. The main differ- 
ence between electrons and ions is that as soon as the sur- 
face collected some electrons, because of the faster bom- 
bardment with electrons than with ions, the surface po- 
tential for ions is the attractive Coulomb potential (most 
probably screened but thats for the following irrelevant). 
Hence, ion surface states develop in the tail of the long- 
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ranged Coulomb potential and thus deep in the sheath 
of the grain, far away from its surface. The microscopic 
processes driving ion energy relaxation and eventually ion 
sticking and desorption are thus not surface- but plasma- 
bound. 

In the microscopic approach presented below, we fo- 
cus on the physics occurring at most a few nm away from 
the boundary. We will therefore not give here a quantita- 
tive treatment of the physisorption kinetics of ions in the 
long-ranged Coulomb potential. However, when it comes 
to the calculation of the surface charge via phenomenolog- 
ical equations connecting the quantum with the classical 
level, we have to make some assumptions about the ion 
dynamics and kinetics. We will then discuss ions qualita- 
tively. The assumptions made for ions, which are some- 
what in conflict with what other people expect 56,57,58 , 
do however not affect the microscopic calculation of s e and 

The outline of this colloquium is as follows. In the 
next section we describe and put into context the sur- 
face model for the charge of a floating dust particle in 
a plasma we developed in |29j because it motivated the 
physisorption-inspired microscopic treatment of electrons 
at plasma boundaries discussed in this colloquium. A qual- 
itative description of the ion kinetics in the vicinity of a 
spherical grain is also included in this section. Section 
3 describes a microscopic model for the interaction of 
electrons with plasma boundaries. Specified to a metal- 
lic boundary, it will then be used to calculate the electron 
sticking coefficient s e and the electron desorption time r e . 
Key issues of the microscopic description of the electron- 
wall interaction (surface potential, coupling to elementary 
excitations of the solid, etc.) will be identified and numer- 
ical results will be presented and discussed. A critique 
of our assumptions is given at the end of section 3 and 
should be understood as a list of to-do's. We close the 
presentation in section 4 with a few concluding remarks. 
Mathematical details interrupting the presentation which 
is meant to be read in order because it successively con- 
structs a case are relegated to three appendices. 

2 Charge of a dust particle in a plasma 

The physisorption-inspired treatment of surface charges 
originated from our attempt to calculate the charge of 
a spherical /xm-sized floating dust particle in a quiescent 
plasma, taking not only plasma-induced but also surface- 
induced processes into account [29] . Here we have to clearly 
distinguish between the assumptions made to construct a 
constituting equation for the surface charge, which by ne- 
cessity has to connect the quantum mechanics occurring 
at the surface with the classical physics determining the 
plasma fluxes, and the assumptions to obtain estimates 
for the surface parameters appearing in this equation. The 
microscopic calculation of the electron surface parameters 
s e and r e presented in the next sections is of course inde- 
pendent of the assumptions about the ion dynamics and 
kinetics as well as the phenomenological nature of the con- 
stituting equation for the surface charge. 



2.1 Rate equations 

First, we will discuss the surface model proposed in [29] 
from the perspective of the rate equations corresponding 
to the elementary processes shown in Fig. [T] Thereby we 
also identify the assumptions, in particular, with respect 
to the surface properties, which are usually made in stan- 
dard calculations of surface charges. 

To be specific let us consider a spherical dust particle 
with radius R. The quasi-stationary charge of the grain is 
given by (we measure charge in units of — e) 

Z p = 4ttR 2 [a e - a,] , (1) 

with electron and ion surface densities, cr e ^, satisfying the 
quasi-stationary (da e ^/dt — 0) rate equations [28] , 

n „ -plasma 1„. ~ ^ „ /o\ 

r, -plasma —1 / \ 

= SiJl - r,- Oi - a R a e Ui , (3) 

where Je 1 I asma , s e ,i, T e ^, and a R denote, respectively, the 
fluxes of electrons and ions hitting the grain surface from 
the plasma, the electron and ion sticking coefficients, the 
electron and ion desorption times, and the electron-ion 
recombination coefficient. 

In order to derive the standard criterion invoked to 
determine the quasi-stationary grain charge, we now as- 
sume, in contrast to what we do in our model [29] (see 
also below), that both electrons and ions reach the sur- 
face of the grain. In that case, both Eq. © and Eq. ^ 
should be interpreted as flux balances on the grain sur- 
face. At quasi-stationarity, the grain is charged to the 
floating potential U. In energy units, U = Z p e 2 /R = 
2Z p Roclb/ R with Rq the Rydberg energy and as the 
Bohr radius. Because the grain temperature ksT s <C U 
the ion desorption rate t^ 1 w 0. Equation ([3]) reduces 

therefore to a R a e Ui = Sij l plasma which transforms Eq. ([2]) 
into s e J| ,lasma = s i jf lasma + T~ 1 a provided a w a e which is 
usually the case. In the standard approach the grain sur- 
face is moreover assumed to be a perfect absorber for both 
electrons and ions. Thus, s e — Sj = 1 and r" 1 = t^ 1 = 0. 
The quasi-stationary charge Z p of the grain is then ob- 
tained from the condition 

j» ias ™(Z p ) = jf las ™(Z p ) , (4) 

where we explicitly indicated the dependence of the plasma 
fluxes on the grain charge. 

Calculations of the grain charge differ primarily in the 
approximations made for the plasma fluxes jg lasma . For 
the repelled species, usually collisionless electrons, the flux 
can be obtained from Poisson's equation and the collision- 
less Boltzmann equation, using trajectory tracing tech- 
niques based on Liouville's theorem and energy and mo- 
mentum conservation 59, 60 , 61 . The flux for the attracted 

1 The rate equations connecting the plasma fluxes j^ 1 f sma 
and surface densities cr e) , with the surface parameters s e ,i, 
T e ,i, and an are phenomenological. They should be derived 
from Boltzmann equations containing surface scattering inte- 
grals which encapsulate the quantum mechanics responsible for 
sticking, desorption, and recombination. 
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largest stable orbit (effective surface for ions) 



Fig. 2. Illustration of the surface model for the charging of 
a dust particle with radius R in a gas discharge. At quasi- 
stationarity, surface charges a e< i bound at r e w R and rt > 
r e , respectively, balance the collection flux s e ^j^ srn ^ with the 
respective desorption flux T~l<r e ,i, where s e ,i and r ej i denote, 
respectively, sticking coefficients and desorption times 29 . 

species, usually collisional ions, is much harder to obtain. 
Unlike the electron flux, the ion flux depends not only on 
the field of the macroscopic body but also on scattering 
processes due to the surrounding plasma, which through- 
out we assume to be quiescent. For weak ion collisionalities 
the charge-exchange enhanced ion flux model proposed 
by Lampe and coworkers [56 ,57,58] is usually used. Its 
validity has been however questioned by Tskhakaya and 
coworkers [62l[63] . We come back to Lampe and coworkers 
approach below when we discuss representative results for 
our surface model. 

Hence, irrespective of the approximations made for the 
plasma fluxes, the standard approach of calculating sur- 
face charges is based on three assumptions about the sur- 
face physics: 

• Both ions and electrons reach the surface, even on the 
microscopic scale. 

• s e = Si = 1 or at least s e = s;. 

— 1 n i i , _1 -plasma 

• r e — or at least r e a e <C Sjjf = ctRO e Ui. 

We basically challenge all three assumptions. 

First, electrons and ions should be bound in surface 
states. Because of differences in the potential energy, mass, 
and size the spatial extension of the electron and ion bound 
states, and thus the average distance of electrons and ions 
from the boundary, is expected to be different. On the mi- 
croscopic scale, electrons and ions trapped to the surface 
should be spatially separated. 

Second, s e = Si is quite unlikely. Usually heavy par- 
ticles, such as ions, couple rather strongly to vibrational 
excitations of the boundary [36,39 . They can thus dissi- 
pate energy very efficiently which usually leads to a large 
sticking coefficient. Light particles, like electrons, on the 
other hand, couple only very weakly to vibrations of the 
solid. On this basis, we would expect s e <§; Si. To what 
extend the coupling to other elementary excitations of the 
boundary (plasmons, electron-hole pairs, ...) can compen- 
sate for the inefficient coupling to lattice vibrations is part 
of our investigations. 



Third, if ions and electrons are indeed spatially sepa- 
rated, the two rate equations should be in fact interpreted 
as flux balances on two different effective surfaces (viz: 
the two closed circles in Fig. [2]). In that case, a_R<7ier e <C 
&e,i/i~ e ,i and the surface charge Z p would be determined 
by balancing on the grain surface the electron desorption 
flux, r" 1 ^, with the electron collection flux, s e jP lasma . 
The corresponding balance of ion fluxes, to be taken on an 
effective surface surrounding the grain, would then yield 
a partial screening charge Zi. Within this scenario, we 
would thus obtain 

Z p = Anrl ■ (sr) e ■ jf^{Z p ) , (5) 
Z 4 = 4^ 2 .( S r) 2 . A plasma , (6) 

with ?' e w R and > r e . 

The surface physics is now encoded in (st)^. These 
products depend on the material and the plasma. They 
could be used as adjustable parameters. A justification 
of the assumptions, however, made in deriving Eqs. (J5J) 
and ((6]) can only come from a microscopic calculation of 
(sr)e,i- 

For electrons, various aspects of this calculation will 
be discussed in the following sections. 



2.2 Semi-microscopic approach 

Before we discuss the complete microscopic calculation of 
s e and T e we summarize the semi-microscopic approach 
taken in Ref . [29] . This prepares the grounds for a micro- 
scopic thinking and demonstrates that Eqs. ((5]) and (J6j) 
give results which compare favorable with experimental 
data. 

The approach we adopted in Ref. [SS] is based on a 
quantum mechanical investigation of the bound states of 
a negatively charged particle in a gas discharge. For that 
purpose, we considered the classical interaction between 
an electron (ion) with charge —e (+e) and a spherical 
particle with radius R, dielectric constant e, and charge 
Z p . The interaction potential contains then a short-ranged 
polarization-induced part arising from the electric bound- 
ary conditions at the grain surface - the classical image 
potential - and a long-ranged Coulomb tail due to the 
particle's charge [641165] . 

The polarization-induced part of the potential will be 
discussed from a quantum- mechanical point of view in ap- 
pendix A. Concerning the Coulomb tail we may add that it 
arises from the interaction between the approaching elec- 
tron and the electrons already residing on the grain. From 
many-body theory it is known that this interaction can 
be rather involved because the attached electrons may re- 
spond dynamically [BH). We neglect this possibility. The 
Coulomb part is then simply the potential of a sphere 
(plane) with charge Z p . This is equivalent to a meanfield 
approximation for the electron-electron interaction. 

Measuring distances from the grain surface in units 
of R and energies in units of U , the interaction energy 
at x = r / R — 1 > Xb, where Xb is a lower cut-off, below 
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which the grain boundary cannot be described as a perfect 
surface anymore, reads 



1 



Ve ' l{x) ~ ± l+x x(l + x) 2 (2 + x) 

1 — £/2x electron 
-1/(1 + x) ion 



(7) 



with £= (e- l)/2(e + l)Z p . 

The second line in Eq. ([7|) is an approximation which 
describes the relevant parts of the potential very well and 
permits an analytical calculation of the surface states. In 
Fig. [3] we plot V et i(x) for a melamine-formaldehyde (MF) 
particle (e = 8, R = 1 fxm, and Z p = 1500) embed- 
ded in a lOOPa neon discharge with plasma density n e = 
rii = 0.39 x 10 9 cto~ 3 , ion temperature ksTi = 0.026 eV, 
and electron temperature ksT e = 6.3 eV [15] . From the 
electron energy distribution, f e (E), we see that the dis- 
charge contains enough electrons which can overcome the 
Coulomb barrier of the dust particle. These electrons may 
get bound in the polarization-induced short-range part of 
the potential, well described by the approximate expres- 
sion, provided they can get rid of their kinetic energy. 
Ions, on the other hand, being cold (see fi(E) in Fig. [3} 



and having a finite radius rf lze /R 



> 10" 



not explore the potential at short distances. For them, the 
long-range Coulomb tail is most relevant, which is again 
well described by the approximate expression. 

Writing for the electron eigenvalue e e = 1 — a e £/4fc 2 
with a e = (e — l)R/4(e + 1)ob and for the ion eigen- 
value e 1 — —a.i/2k 2 with ai — miRZ p /m e aB, where m e 
and rrii are the electron and ion mass, respectively, the 
radial Schrodinger equations with the approximate poten- 
tials read 



d 2 u e 



dx 2 



-g+Ve,i(x) 



1(1 + 1) 



(l + x) 2 



(8) 



where V e (x) — 2a e /x and Vi(x) = 2a.i/(\ + x). 

For bound states, the wavefunctions have to vanish for 
x — > oo. The boundary condition at xt depends on the 
potential for x < Xb, that is, on the potential within the 
solid (which is different for electrons and ions). Match- 
ing the solutions for x < Xb and x > Xb at x = Xb 
leads to a secular equation for k. Ignoring the possibil- 
ity that electrons and ions may also enter the solid, we set 
V e ,i(x < Xb) — oo with Xb — for electrons and Xb = xf ze 
for ions. For electrons we thereby restrict ourselves to 
weakly bound polarization-induced surface states, neglect- 
ing strongly bound crystal-induced surface states which, 
in general, may also occur [67j . As explained in the next 
section, we expect them to be of minor importance for 
physisorption of electrons. 

The electron Schrodinger equation with the hard bound- 
ary condition at z — is equivalent to the radial Schrodinger 
equation for the hydrogen atom. Hence k is an integer 
n. Because (for bound electrons) x -C 1 and a e >■ 1, 
the centrifugal term is negligible. We consider therefore 
only states with 1 = 0. The eigenvalues are then e e n = 




Fig. 3. Left panel: Potential energy for an electron (ion) in the 
field of a MF particle (R = 1 /im, Z = 1500) [15] and represen- 
tative probability distributions, |m(x)| 2 , shifted to the binding 
energy and maxima normalized to one. Dashed lines denote the 
potentials used in the Schrodinger equations. Note, the finite 
ion radius rf 2e ~ A forces the ion wavefunctions to vanish at 
x ~ 10~ 4 . Right panel: Bulk energy distribution functions for 
the lOOPa neon discharge hosting the particle [JST: ksT e = 
6.3eV, k B Ti = 0.026eV, and n e = m = 0.39 x 10 9 cm" 3 . 



1 — a e ^/4n 2 and the wavefunctions read 

u nfi( X ) ~ v nfi{z) 

= zexp(-z72)(-)- 1 (n - l)!!^) (9) 

with z = 2a e x/n and Ln\z) associated Laguerre polyno- 
mials. 

The probability densities |u^ (a;)| 2 for the first three 
states are plotted in Fig. [3] As can be seen, electron sur- 
face states are only a few Angstroms away from the grain 
boundary. At these distances, the spatial variation of V e (x) 
is comparable to the de-Broglie wavelength of electrons 
approaching the particle. More specifically, for ksT e = 
6.3 eV, X d e B /R w |K/K'I ~ 10_4 - Hence, the trapping of 
electrons at the surface of the particle has to be described 
quantum- mechanically. 

The solutions of the ion Schrodinger equation are Whit- 



taker functions, u l k ^x) = VFfc,i+i/2( x ) with x = 2tVi(l + 
x)/k and k determined from u\ i(xf lze ) = 0. However, 
since k 3> 1 and x 3> 1, it is very hard to work directly 
with W k ^ + i/ 2 {x). It is easier to use the method of compar- 
ison equations [68J and to construct uniform approxima- 
tions for u\ i (x) with the radial Schrodinger equation for 
the hydrogen atom as a comparison equation. The method 
can be applied for any I. Here we give only the result for 
/ i): 



u k o( x ) ~ v n,o( z )/ V dz/dx 



(10) 



with v n $(z) defined in Eq. ^ and z = 2otiz(x)/n. The 
mappings z(x) and k(n) can be constructed from the phase- 
integrals of the two Schrodinger equations. 

In Fig. [1 we show K i0 (x)| 2 for fc(300) and fc(30000). 
Note, even the fc(30000) state is basically at the bottom 
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of the potential. This is a consequence of a, 3> 1 which 
leads to a continuum of states below the ion ionization 
threshold at e = 0. We also note that |wjL n ) (a;)| 2 peaks 
for n» 1 just below the turning point. Hence, except for 
the lowest states, which we expect to be of little impor- 
tance, ions are essentially trapped in classical orbits deep 
in the sheath of the grain. This will be also the case for 
/ > 0. That ions behave classically is not unexpected be- 
cause for UbTi — 0.026 eV their de-Brogfie wavelength is 
much smaller then the scale on which the potential varies 
for x > 10~ 3 : Xf B /R « 10~ 5 < \Vi/V?\ » 1. Thus, the 
interaction between ions and the particle is classical. 

Nevertheless it can be advantageous to describe ions 
quantum-mechanically and to use the method of compar- 
ison equations, which is an asymptotic technique, to per- 
form the calculation in the semiclassical regime. Since the 
ion dynamics and kinetics is beyond the scope of this pa- 
per, we do not give more mathematical details about the 
solution of the ion Schrodinger equation. We mention how- 
ever that many years ago Liu [69] pursued a quantum- 
mechanical description of the collisionless ion dynamics 
around electric probes. But he found no followers. 

A model for the charge of the grain which takes sur- 
face states into account can now be constructed as fol- 
lows. Within the sheath of the particle, the density of 
free electrons (ions) is much smaller than the density of 
bound electrons (ions). In that region, the quasi-stationary 
charge (again in units of — e) is thus approximately given 

by 



Z[x) = 4vri? 3 / dx'(l 



(11) 



with x < \f = y/kBTi/4:TTe 2 n.i, the ion Debye length, 
which we take as an upper cut-off, and n b e i the density of 
bound electrons and ions. For the plasma parameters used 
in Fig. [3j Af « 60/im. The results for the surface states 
presented above suggest to express the density of bound 
electrons by an electron surface density: 



n b e {x) « a e 8(x — x e )/R 



(12) 



with x e xi, ~ and a e the quasi-stationary solution of 
of Eq. ([2]) without the recombination term. Equation ([2]) 
is thus still interpreted as a rate equation on the grain sur- 
face. We will argue below that once the grain has collected 
some negative charge, not necessarily the quasi-stationary 
one, there is a critical ion orbit at Xi ~ 1 — 10 ^> x e which 
prevents ions from hitting the particle surface. Thus, the 
particle charge obtained from Eq. (jlip is simply Z p = 
Z(x e < x < x^. Inserting Eq. (fl~2"j) into Eq. (fTTjl and in- 
tegrating up to x with x e < x < Xi leads to Eq. ([5]), the 
expression for the particle charge deduced from the rate 
equations ^ and ([3]) under the assumption that ions do 
not reach the grain surface on the microscopic scale. 

For an electron to get stuck at (to desorb from) a 
surface it has to loose (gain) energy at (from) the sur- 
face [361 . This can only occur through inelastic scattering 
with the grain surface. To calculate the product (sr) e re- 
quires therefore a microscopic description of energy relax- 
ation at the grain surface. This will be discussed in the 



next section. In Ref. [29] we invoked the phenomenology 
of reaction rate theory and approximated (sr) e by 



(sr)e 



k B T t 



exp 



(13) 



where h is Planck's constant, T s is the surface tempera- 
ture, and -Eg is the electron desorption energy, that is, the 
binding energy of the surface state from which desorption 
most likely occurs [361 . The great virtue of this equation is 
that it relates a combination of kinetic coefficients, which 
depend on the details of the inelastic (dynamic) interac- 
tion, to an energy, which can be deduced from the static 
interaction alone. Kinetic considerations are thus reduced 
to a minimum. They are only required to identify the rel- 
evant temperature and the state from which desorption 
most probably occurs. In the next section we will show, 
for a particular model, how Eq. (|13p can be obtained from 
a microscopic theory. Its range of validity will then become 
also clear. 

Equation ([5]) is a self-consistency equation for Z p . Com- 
bined with Eq. (| 13(1 , and approximating the electron flux 
^-plasma f rom pi asma by the orbital motion limited flux, 



3 



OML 



= n e ^/kj^I\J^jTm e eyi])[—Z p e 2 /RkBT e ] , (14) 



which is reasonable, because, on the plasma scale, elec- 
trons are repelled from the grain surface, the grain charge 
is given by 



A-kR' 



k B T s 



T 'j? ML (z p ) . 



(15) 



Thus, in addition to the plasma parameters n e and T e , the 
charge depends on the surface parameters T s and E^. 

Without a microscopic theory for the inelastic electron- 
grain interaction, a plausible estimate for Ef has to be 
found from physical considerations alone. Since by neces- 
sity the electron comes very close to the grain surface (see 
Fig. [3]) it will strongly couple to elementary excitations of 
the grain. Depending on the material these may be bulk 
or surface phonons, bulk or surface plasmons, or internal 
electron- hole pairs. For any realistic description of the po- 
tential for x < xj, the electron wavefunction leaks into the 
solid, the electron will therefore quickly relax to the lowest 
surface bound state. The microscopic model for electron 
energy relaxation at metallic boundaries presented in the 
next section turns out to even work for an infinitely high 
barrier. Taking the n — 1 state for V e ^(x < Xb) — oo as 
an approximation to the lowest surface bound state, it is 
reasonable to expect 



(l-e\)U = 



Rnfe-1 



16 



1 



(16) 



which, for an MF particle with e = 8, leads to E^ 0.5eV. 
The particle temperature cannot be determined in a sim- 
ple way. It depends on the balance of heating and cooling 
fluxes to-and-fro the particle and thus on additional sur- 
face parameters [70] ■ We use T s therefore as an adjustable 
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?>.5 1 1.5 

Radius [10" 6 m] 

Fig. 4. Radius dependence of the charge of a MF particle in 
the bulk of a neon discharge at p = lOOPa [15| . The plasma 
parameters are the same as in Fig. [3] The solid line denotes 
the charges deduced from Eq. (I15|l and the dashed line gives 
the charges obtained from j° ML = j° ML + jf x with <t cx = 
10- 14 cm" 2 . 

parameter. To reproduce, for instance, with Eq. (|15p the 
charge of the particle in Fig. [3] T s = 370 K implying 
(sr) e PS 10" 6 s. 

In Fig. 3] we plot the radius dependence of the charge 
of a MF particle in the lOOPa neon discharge specified in 
the caption of Fig. [3] More results are given in [29]. Since 
the plasma parameters are known the only adjustable pa- 
rameter is the surface temperature. Using T s = 370K we 
find excellent agreement between theory and experiment. 
For comparison we also show the charges obtained from 
Eq. Q , approximating the ion plasma flux by 

i pla, ma = i OML +i cX j 

where 

jP ml = WM^vrm^l + Z p e 2 /Rk B T t ] (18) 
is the orbital motion limited ion flux and [15] 

3 T = ^(OTAf /l CK )^k B T l /2jrm 1 (Z p e 2 /Rk B Ti) 2 (19) 

is the ion flux originating from the release of trapped 
ions due to charge-exchange scattering as suggested by 
Lampe and coworkers [56l[57j|58] . The scattering length 
'cx = (fcx^g) 1 with <t cx = 10~ 14 cm 2 the scattering cross 
section and n g — p/k B T g the gas density. Clearly, the ra- 
dius dependence of the grain charge seems to be closer to 
the nonlinear dependence obtained from Eq. (fT5|) than to 
the linear dependence resulting from 

j? ML = j? ML +jr, (20) 

indicating that the surface model we propose captures at 
least some of the physics correctly which is responsible for 
the formation of surface charges. 



In order to derive Eq. (fT5| from Eq. (fTTj) we had to as- 
sume that once the particle is negatively charged ions are 
trapped far away from the grain surface. Treating trapping 
of ions in the field of the grain as a physisorption pro- 
cess suggests this assumption, which is perhaps counter- 
intuitive. Similar to an electron, an ion gets bound to the 
grain only when it looses energy. Because of its low en- 
ergy and the long-range attractive ion-grain interaction, 
the ion will be initially bound very close to the ion ion- 
ization threshold (see Fig. [3]). The coupling to the ele- 
mentary excitations of the grain is thus negligible and 
only inelastic processes due to the plasma are able to 
push ions to lower bound states. Since the interaction is 
classical, inelastic collisions, for instance, charge-exchange 
scattering between ions and atoms, act like a random 
force. Ion energy relaxation can be thus envisaged as a 
de-stabilization of orbits. This is in accordance to what 
Lampe and coworkers assume [56,57,58 . In contrast to 
them, however, we |29j expect orbits whose spatial ex- 
tension is smaller than the scattering length to be sta- 
ble because the collision probability during one revolution 
becomes vanishingly small. For a circular orbit, a rough 
estimate for the critical radius is 

n = R(l + x t ) = (2 7 rcr cx n 9 )- 1 (21) 

which leads to Xi ~ 5.7 3> x e ~ when we use the param- 
eters of the neon discharge of Fig.[3]and cr cx = 10 -14 cm 2 . 

Although the approach of Lampe et al. [Mi pTJoSj shows 
a pile-up of trapped ions in a shell of a few [im radius 
enclosing the grain, they would not expect a relaxation 
bottleneck. This point can be only clarified with a de- 
tailed investigation of the ion dynamics and kinetics in 
the vicinity of the grain, including electron-ion recombi- 
nation. As mentioned before, despite the classical charac- 
ter of the ion dynamics, a quantum-mechanical treatment, 
similar to the one we will present in the following sections 
for electrons, is possible and perhaps even advantageous 
because it treats closed (bound surface states) and open 
ion orbits (extended surface states) on the same footing. 
In addition, energy barriers due to the angular motion 
are easier to handle in a quantum- mechanical context. In 
fact, Lampe and coworkers neglect these energy barriers 
whereas Tskhakaya and coworkers [62 , 63J believe that this 
approximation overestimates j 4 cx . In reality, they claim, 

is much smaller. If this is indeed the case, the condi- 
tion ML = jf ML + jf x would yield charges which are 
much closer to the orbital-motion limited ones and thus 
far away from the experimentally measured charges. 

Pushing the assumption of a critical ion orbit even 
further, we assumed in [29] that all trapped ions can be 
subsumed into a single effective orbit as shown in Fig. [2] 
We then obtained an intuitive expression for the number 
of ions accumulating in the vicinity of the grain, that is, 
for its partial screening charge. For that purpose we mod- 
elled the ion density n\ accumulating in the vicinity of 
the critical orbit by a surface density <Ji which balances 
at Xi the ion collection flux Sijf lasma with the ion des- 
orption flux r" 1 ^,;. Mathematically, this gives rise to a 
rate equation similar to ([3]) , with the recombination term 
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neglected and interpreted as a rate equation at r = r,. 
Although Eq. (fl3|) assumes excitations of the grain to be 
responsible for sticking and desorption we expect a simi- 
lar expression (with E%, T s replaced by Ef , T g ) to control 
the density of trapped ions. Integrating (jlip up to x with 
Xi < x < Xf we then obtain Z(xi < x < Xf) = Z p — Z% 
with 

Zi = 4t:R 2 (1 + xtf — e&W'B^j? (22) 
k B T g 

the number of trapped ions. Since the critical orbit is near 
the sheath-plasma boundary, it is fed by the Bohm ion flux 

if = Ofimy/kBTe/mi . (23) 

The ion desorption energy is the negative of the binding 
energy of the critical orbit, 

Ef(Z p ) = -Vi(xi)U(Z p ) = 4Tr<j cx a B n g Z p Ro , (24) 

and depends strongly on Z p and Xi. For the situation 
shown in Fig. [3l we obtain Ef « 0.39eT^ and (sr)i w 
10~ 8 s when we use T g = T s = 370 K, the particle tem- 
perature which reproduces Z p w 1500. The ion screen- 
ing charge is then Zi ~ 12 <C Z p which is the order 
of magnitude expected from molecular dynamics simula- 
tions [71] . Thus, even when the particle charge is defined 
by Z(xi < x < Xf ) it is basically given by Z p . 

From the surface model we would expect (sr) e ~ 10~ 6 s 
to produce particle charges Z p of the correct order of mag- 
nitude. Since the particle temperature T s is unknown, it 
can be used as an adjustable parameter. The calculated 
Z p can thus be always made to coincide with the mea- 
sured charge. The particle temperature has to be of course 
within physically meaningful bonds. Recently, the particle 
temperature (but unfortunately not the particle charge) 
has been measured |72j . There is thus some hope that in 
the near future Z p and T s will be simultaneously mea- 
sured. Finally, let us point out that, because ions are in 
our model bound a few microns away from the surface, we 
obtain (sr)j < {sr) e , in agreement with the phenomeno- 
logical fit performed in [28 . 



3 Physisorption of electrons 

In the previous section we described a microscopic, physi- 
sorption-inspired model for the charging of a dust parti- 
cle in a plasma which avoids the unrealistic treatment of 
the grain as a perfect absorber. Within this model the 
charge and partial screening of a dust particle can be cal- 
culated without relying on the condition that the total 
electron plasma flux balances on the grain surface the 
total ion plasma flux. Instead, two flux balance condi- 
tions are individually enforced on the two effective sur- 
faces shown in Fig. [2] (solid circles). The quasi-stationary 
particle charge Z p is then given by the number of elec- 
trons "quasi-bound" in the polarization potential of the 
grain and the screening charge Zi is approximately given 



by the number of ions "quasi-trapped" in the largest sta- 
ble closed ion orbit (which defines an effective surface for 
ions and subsumes, within our model, all trapped ions into 
a single effective orbit). 

The physisorption kinetics at the grain boundary, that 
is, the sticking in and the desorption from (external) sur- 
face states due to inelastic scattering processes, is en- 
coded in the products (sr) e ^ which we approximated by 
phenomenological expressions of the form (|13p . For elec- 
trons, we now take a closer look at what happens on the 
surface microscopically. First, we will discuss the micro- 
physics qualitatively. Then we will perform an exploratory 
quantum-mechanical calculation of s e and r e using a sim- 
ple one-dimensional model for the electronic properties of 
the surface which allows us to do large portions of the cal- 
culation analytically. Finally, we will critically assess the 
results of the calculation turning thereby its shortcomings 
into a list of to-do's. 

In principle, trapping and de-trapping of ions in the 
surface-induced Coulomb potential of the grain can be 
also understood as a physisorption process. However, the 
quantum-mechanical approach we will use for electrons 
has then to be pushed to the semi-classical regime ap- 
propriate for ions. In addition, not surface- but plasma- 
based inelastic scattering processes will turn out to control 
ion energy relaxation. Although conceptually very close, 
mathematically the calculation of Si and t% is quite dif- 
ferent from the calculation of s e and r e . It is therefore 
beyond the scope of this paper. In the concluding section 
we may however add a few remarks about ions. 



3.1 Qualitative considerations 

The surface of a /Ltm-sized grain is large enough to contain 
sizeable spatial regions (facets) isomorphous to crystallo- 
graphic planes. Except specific features arising from the 
finite extend of the facets, whose influence diminishes with 
the facet size, the electronic properties of the facets resem- 
ble the electronic properties of (infinitely extended) planar 
surfaces. In particular, like ordinary surfaces, facets should 
support surface states to which electrons approaching the 
grain from the plasma may get bound and then re-emitted 
when they dynamically interact with the elementary exci- 
tations of the grain. 

Each facet may give rise to two types of surface states: 
(i) Crystal-induced surface states due to the abrupt ap- 
pearance (from the plasma electron's point of view) of a 
periodic potential inside the grain and (ii) polarization- 
induced image states, on which the considerations of the 
previous section were based. Compared to the binding en- 
ergy of image states, the binding energy of crystal-induced 
surface states is very large. Instead of a few tenth of an 
electron volt, it is typically a few electron volts. As a re- 
sult, the center of gravity of crystal- induced surface states 
is much closer to the surface than the center of gravity of 
image states. 

Based on the experimental results of [UJ we show in 
Fig. [5] as an example the schematic electronic structure 
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Cu(100) Cu(llO) Cu(lll) 




Fig. 5. Schematic electronic structure after |47| for a copper 
(100), a copper (110), and a copper (111) surface, respectively, 
for the lateral momentum where the projected energy gap is 
largest: \Q\ = for Cu (100) and Cu (111); \Q\ = 1.2A" 1 
in [110] direction for Cu (110). Shaded areas denote the pro- 
jected bulk band structure and thick solid lines indicate crystal- 
induced (n = 0) and polarization-induced (n > 1) surface 
states. The zero of the energy axis is the Fermi energy and 
the vacuum level is given by the dashed line. 

of three copper surfaces, respectively, for the lateral mo- 
mentum where the projected energy gap is largest. The 
electronic structure for a given orientation changes with 
momentum (not shown) but for all orientations, and that 
is the point we want to make, surface states exisiQ, in ad- 
dition to projected bulk states, and may thus participate 
in a physisorption process. For dielectric surfaces the elec- 
tronic structure is quite similar although the details and 
physical origin of the states is different [67] . 

An ab-initio modeling of surface states is complex and 
computationally expensive, even for planar surfaces (see 
for instance [73]). Fortunately, the essential physics can be 
understood within simple one dimensional models which 
assume the potential energy to vary only normally to the 
surface (z— direction) as illustrated in Fig. [6] Inside the 
material (z < 0) the potential has the periodicity of the 
crystal. It may thus lead to an energy gap on the sur- 
face. Outside the material, the potential gives rise to a 
barrier which merges at large distances with the asymp- 
totics of the image potential V p (z) ~ —1/z. Its physical 
origin are exchange and correlation effects which, on the 
one hand, contribute to the confinement of electrons inside 
the material and, on the other hand, cause the attraction 
of external electrons to the surface. A simple microscopic 
model for the image potential 41,42] is given in appendix 
A. 

The situation shown in Fig. [5] is the most favorable 
one for physisorption of electrons. The vacuum (plasma) 
potential, which is the zero of the energy scale, is in the 
middle of a large energy gap. Four main classes of states 

2 More precisely, surface states exist in some parts of the 
surface Brillouin zone [47] . 



can then be distinguished: (i) Volume states periodic in- 
side the material and exponentially decaying into the vac- 
uum (plasma). They exist for energies where bulk states 
are also allowed. Close to band edges they may have an in- 
creased weight near the surface in which case they are sur- 
face resonances, (ii) Bound surface states, that is, states 
decaying exponentially into the material and the vacuum. 
They appear in regions of negative energies where bulk 
states are absent: Weakly bound image states close to the 
vacuum potential and strongly bound crystal-induced sur- 
face states close to the Fermi energy. Crystal-induced sur- 
face states may have tails on the material side strongly 
oscillating with the crystal periodicity, while the tails of 
image states may only weakly respond to the crystal po- 
tential, (iii) Unbound surface states for positive energies 
inside the gap. They are free on the vacuum and bound on 
the material side. The periodic crystal potential may also 
not affect these states very much, (iv) States which are 
free on both sides. Inside the material they oscillate with 
the lattice periodicity while outside the material their os- 
cillations have to fit the surface potential. In the vicinity 
of the surface this class of states may also have a peak. 

Of particular importance for sticking and desorption 
are transitions between bound and unbound surface states 
due to inelastic scattering with elementary excitations of 
the boundary. The elementary excitations can be phonons, 
plasmons, and electron-hole pairs. The latter two cases are 
excitations involving volume states. 

The potential plotted in Fig. [6] is for an uncharged 
surface. An electron approaching a plasma boundary is of 
course also subject to the Coulomb repulsion due to the 
electrons already residing on the surface. In the mean- 
field approximation, however, this repulsion leads only to 
a barrier whose height is the floating energy U. Only an 
electron with an energy larger than U has a chance to 
come close enough to the surface to feel the attractive 
part of the potential. For an electron bound in this part 
of the potential, on the other hand, the Coulomb barrier 
merely sets the ionization threshold. Thus, as long as the 
Coulomb repulsion is treated in meanfield approximation, 
the Coulomb term drops out from the considerations pro- 
vided we shift the zero of the energy axis to U, that is, 
by simply measuring energies with respect to the floating 
energy (Coulomb barrier) of the surface. If U falls inside 
an energy gap of the boundary the situation is similar to 
the one depicted in Fig. [6] 

3.2 Simplified planar microscopic model 

Ideally, a microscopic calculation of s e and r e for a spher- 
ical grain would be based on a three-dimensional first- 
principle electronic structure of the grain surface. 

In view of the discussion of the previous subsection an 
estimate for the grain's s e and r e may be however also 
obtained by the following strategy which is most proba- 
bly simpler because it allows to incorporate existing (one- 
dimensional) empirical pseudo-potentials for planar sur- 
faces: (i) Identify the facets on the grain surface and ne- 
glect, in a first approximation, the finite lateral exten- 
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z=0 



vacuum (plasma) 



Fig. 6. Schematic drawing of the potential energy at a surface 
(plasma boundary) such as copper (100) or copper (110) and 
representative wavefunctions for volume states (i), bound sur- 
face states (ii), unbound surface states (iii), and free states (iv). 
Shaded areas denote projected bulk states, Ef is the Fermi en- 
ergy of the solid, E — is the vacuum (i.e. floating) level, and 
fe{E) is the plasma electron's Boltzmann distribution function. 
The dashed lines indicate the approximate potentials defin- 
ing the simplified planar model of subsection 13.21 on which the 
calculation of s e and r e is based we describe, respectively, in 
subsection 13.31 and 13.41 



sion of the facets, that is, work with plane waves or Bloch 
functions in the lateral dimensions, (ii) Use empirical one- 
dimensional potentials for planar surfaces [53] to calcu- 
late for each facet separately bound and unbound sur- 
face states, (iii) Identify the channels for electron energy 
relaxation and set up, again for each facet separately, 
a quantum-kinetic scheme for the calculation of s e and 
r e . (iv) Use an appropriate macroscopic spatial averaging 
scheme to obtain an estimate for the grain's s e and r e . 

Despite its approximate nature this strategy is still 
demanding. To work it out for a realistic grain is surely 
beyond the scope of this colloquium. In the exploratory 
calculation of s e and r e presented below we focused there- 
fore on a single, infinitely extended facet, that is, on a 
planar surface, whose electronic structure we moreover did 
not deduce from an empirical pseudo-potential but from a 
model potential which is amenable to analytical treatment 
while at the same time it retains the essential physics. 

Quite generally, the probability with which an electron 
approaching from the plasma halfspace z > the plasma 
boundary at z — ends up in a bound surface state (stick- 
ing) , or with which an electron bound to the surface ends 
up in a free state (desorption) can be obtained from a 
Hamiltonian, 



Table 1. Dielectric constant e, Debye Energy fcsTb, and the 
energy separations between the four lowest image states for 
different materials. 



material 


e 


k B T D [eV] 


AE 21 [eV] 


AE 43 [eV] 


Cu 


00 


0.03 


0.64 


0.12 


Si 


12 


0.057 


0.46 


0.09 


graphite 


12 


0.19 


0.46 


0.09 


Ceo 


4.5 


0.016 


0.26 


0.05 



where the first term describes the electron motion in the 
static surface potential, the second term denotes the free 
motion of the elementary excitations of the boundary con- 
trolling electron energy relaxation at the boundary and 
thus physisorption of electrons, and the third term is the 
coupling between the two. 

It is advantageous to express the Hamiltonian (f2"S"]) in 
terms of creation and annihilation operators for the (ex- 
ternal) electron as well as the (internal) elementary exci- 
tations. For that purpose we use the basis in which H e is 
diagonal, that is, the eigenstates of the static surface po- 
tential V(z).0 Writing r = (R, z) for the electron position, 
the Schrodinger equation defining these states reads 



2m„ 



A + V(z) W Qq {R,z) =E Qq V Qq (R,z) . (26) 



The lateral motion is free and can be separated from the 
vertical one. Hence, 

& Qq (R, z ) = ^7j G MiQ ■ R]%{z) , (27) 



with A the area of the surface, which is eventually made in- 
finitely large, Q = (Q x , Qy) a two-dimensional wavevector 
characterizing the lateral motion of the electron and ipq( z ) 
the wavefunction for the vertical motion which satisfies the 
one-dimensional Schrodinger equation (viz: Eq. ([5])) 



dz~^ z) 



2m e 



E q - V{z) 



(28) 



with E q = Eq q — h Q /2m e . The quantum number q is 
an integer n for bound and a wavenumber k for unbound 
surface states. In this basis, 



H e - E Ql C Qq C Ql 

Qq 



(29) 



where Cq creates an electron in the surface state WQ q 
with energy 



E Qq = h 2 Q 2 /2m e + E q 



(30) 



The second and third term in (|25p depend on the kind 
of elementary excitations responsible for energy relaxation 



H = H e + H s +H e 



(25) 



3 In a realistic calculation V(z) should be an empirical 
pseudo-potential of the type proposed in [53] . 
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and hence on the material. For dielectric materials, such 
as graphite or silicon, the coupling to vibrational modes 
is most probably the main driving force for physisorption 
of electrons. In particular, lattice vibrations should play 
an important role. Their energy scale is the Debye energy 
ksTu- For most dielectrics ksTo is smaller than the en- 
ergy spacing of the lowest surface states. For image states 
typical energy separations are given in table [TJ When 
crystal-induced surface states or dangling bonds [67]) are 
also included the situation does not change much, it may 
be even worse. Multiphonon processes could thus signif- 
icantly affect physisorption of electrons at dielectric sur- 
faces making it a very interesting problem to study. 

For metals, on the other hand, electronic excitations, 
most notably electron- hole pairs, provide an efficient chan- 
nel for electron energy relaxation [3"T1|40| . They are not 
created across a large energy gap, as in dielectrics, where 
they are therefore unimportant, but with respect to the 
Fermi energy of a partially filled band. In metals electron- 
hole pairs can be excited even at room temperature. Ph- 
ysisorption of electrons at metallic plasma boundaries, 
whose temperatures are typically not much higher than 
room temperature, is thus most likely controlled by the 
coupling to electron-hole pairs. 

The Fermi energy of a metal is inside a band. Electron- 
hole pairs are thus excitations involving volume states. Ig- 
noring, in a first approximation, exchange between these 
states and (bound and unbound) surface states, which 
should be small because the states are spatially separated 
(see Fig.[BI), electrons occupying these two classes of states 
can be approximately treated as two separate species: Ex- 
ternal and internal electrons, where the latter are respon- 
sible for energy relaxation of the former. 

Specifically for a metallic plasma boundary, and we 
will restrict the calculation of s e and r e presented in the 
next two subsections to this particular case, H s is thus the 
Hamiltonian of a non-interacting gas of electronic quasi- 
particles with Fermi energy Ep. Hence, 



Kk 



KkD Kk DKk 



(31) 



with D Kk creating an internal electron in a quasi-particle 
state 



$Kk(R, «) = 4t e MiK ■ R]Mz) (32) 



with energy 



E 



Kk 



h 2 K 2 

2m„ 



Eu 



(33) 



In the above expressions we ignored the periodic crystal 
potential inside the material. It could be taken into ac- 
count using for the lateral motion Bloch states instead of 
plane waves and effective electron masses, possibly differ- 
ent for the lateral and vertical motions, instead of the bare 
electron mass. 

The function <pk{z), describing the vertical motion of 
an internal electron, obeys a one-dimensional Schrodinger 



Table 2. Fermi energy Ef, Fermi wavenumber and screen- 
ing wavenumber (fc s )buik for various metals |74] , 



metal 


E F [eV] 


kplA- 1 ] 


(fcs)bulk/fcF 


Ag 


5.49 


1.20 


1.42 


Cu 


7.0 


1.36 


1.33 


Al 


11.7 


1.75 


1.17 



equation: 



d 2 



2m e 



E k - V(z) 



(34) 



Strictly speaking, the potential V(z) = V(z). But the spa- 
tial parts of the potential determining, respectively, sur- 
face and volume states are different. Working conceptually 
with two separate potentials gives us the flexibility to inde- 
pendently extend the relevant parts of the potential such 
that the calculation of surface and volume states can be 
most easily performed while the essential physics is kept 
(see dashed lines in Fig. [6] and below for the particular 
form of the approximate potentials). 

For a metallic boundary, the interaction part H es of 
the Hamiltonian ([25]) describes the interaction between 
internal and external electrons. Anticipating a statically 
screened Coulomb interaction. 



Hen = 



- E 

QqQ'q' KkK'k 



V^K^ Qg C QW D^ Kk D Klk ,{35) 



with 



V, 



Kk 27re 2 5{Q Q' + K- K') 



Q'q'K'k' 



a 2 VW\ 
i q q UQ-Q') 



(Q - Q') 



(36) 



and 



dzdz'r q (zW k (z')e~ dlz - z ' l 0k'(z')Mz),(37) 



where d = ^fk 2 + Q 2 and k s = (/c s ) sur f ace is the screening 
wavenumber at the surface. Little is known about this 
parameter except that it should be less then the bulk 
screening wavenumber (fc s )buik because the electron den- 
sity in the vicinity of the boundary is certainly smaller 
than in the bulk. In [37] it was for instance argued, based 
on a comparison of experimentally and theoretically ob- 
tained branching ratios for positron trapping at and trans- 
mission through various metallic films that (fc s ) sur f ace ~ 
0.6(fc s )buik- Bulk screening wavenumbers for some metals 
are given in table [2] 

The Hamiltonian (|25[) with H e , H s , and H es respec- 
tively given by ([29]) , (|3l"]) , and ([35]) can be used to calculate 
the transition rate from any initial surface state V'QV to 
any final surface state ^Qq- For the sticking process the 
initial state belongs to the continuum of surface states 



12 



F. X. Bronold et al.: Physisorption kinetics of electrons at plasma boundaries 



K' k' 



Kk 



Qq 



Q'q' 



Fig. 7. Diagrammatic representation of the golden rule (|38|l 
for a transition from the surface state (Q'q') to the surface 
state (Qq) (solid lines) via scattering on an internal electron 
(dashed line) which can be interpreted as the coupling to in- 
ternal electron-hole pairs. The wavy line denotes the screened 
Coulomb interaction between internal and external electrons 
and the box symbolizes the dynamic, that is, inelastic electron- 
metal interaction. 



and the final state is a bound surface state while for the 
desorption process it is vice versa. In lowest order pertur- 
bation theory (see Fig. [Jj , the rate is given by the golden 
rule, 



MQg,QY) = f EEI^V^I 2 

KK' kk' 

x n F (E K , k <)[l - n F (E Kk )\ 

x 6(E QW + E K , k , - E Qq - E Kk ) ,(38) 

where n F (E) = l/(exp[(E - E F )/k B T s ] + 1) is the Fermi 
distribution function for the metal electrons with Fermi 
energy E F and temperature T s . 

To calculate the matrix element ([37)) we need the so- 
lutions of the Schrodinger equations ([2"8)l and Ph- 
ysisorption of electrons involves transitions between bound 
and unbound surface states. The matrix elements for these 
transitions are large when the spatial overlap between the 
initial and final states is large. With unbound surface 
states inside the gap, image states, that is, bound sur- 
face states close to the zero of the energy axis (see Fig- E]), 
have the largest overlap. Crystal-induced surface states, 
having most weight in regions where the weight of un- 
bound surface states is very small, give rise to a smaller 
overlap and are thus less important. We neglect therefore 
crystal-induced surface states and replace V(z) in (12"5)) by 



V(z) 



where 



oo for z < 
V p (z) for z > 



(39) 



(40) 



is the classical image potential. As explained in appendix 
A, V p can be understood in terms of virtual surface plas- 
mon excitations [4Tj|4"2"ll4"3"] . We thus calculated the surface 
states as if the energy gap on the surface were infinitely 
large. The solutions of are then Whittaker functions 
which vanish for z < (see Appendix B) and the required 
matrix elements can be obtained analytically. 



continuum states 




"quantum-mechanical 
boundary layer" 



metal 



Fig. 8. Illustration of the microphysics for a plasma electron 
approaching a metallic boundary. The wavy line indicates the 
surface mode responsible for the attractive polarization poten- 
tial ("image potential") and L is the width of the "boundary 
layer" where a quantum- mechanical calculation applies. The 
electron looses (gains) energy due creation (annihilation) of 
electron-hole pairs in the metal. Due to these processes it may 
get trapped in (escape from) the bound states of the polar- 
ization potential. In other words, it may get stuck at (desorb 
from) the plasma boundary. W and E F are, respectively, the 
work function and the Fermi energy of the metal. 



As far as volume states required for the construction of 
internal electron-hole pairs are concerned, we followed [37J 
|40] and calculated these states as if the work function W 
of the metal where infinite. Measuring moreover energies 
inside the material from the average of the crystal poten- 
tial and neglecting the oscillations of the potential, that 
is, treating the metal boundary as a jellium halfspace, 



V(z) 



for z < 
oo for z > 



(41) 



The wavef unctions <j) k (z) vanish then for z > and are 
standing waves for z < 0. Using box-normalization, 



(42) 



h(z) = \j — siu(kz) 



leading to E k = h 2 k 2 /2m e with k = irn/L and n > 1 
an integer. In the final expressions for s e and r e we use 
L — > oo making k continuous. 

The physical content of the simplified planar model is 
summarized in Fig. [51 It will be used in the next two sub- 
sections to calculate, respectively, s e and r e for a metallic 
plasma boundary. Due to the approximate potentials (|39[) 
and (JU), external and internal single electron wavefunc- 
tions vanish in complementary halfspaces. As a result, the 
matrix element (I37[) factorizes, 



r( 2 ), 



(43) 
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with 



/>00 

I^{Q)= dzeM-^k 2 s +Q 2 W q {z)Tp q >{z) , (44) 
Jo 

roc 

4 2 J'(Q)= dzexp[-z^/k 2 a + Q*]ct>t(-z)4>v(-z) (45) 



to be calculated explicitly in appendix B. 

A rigorous calculation of s e and r e , taking for instance 
into account that sticking and desorption occur on differ- 
ent timescales |35j . should be based on quantum- kinetic 
master equations for the time-dependent occupancies of 
the surface states ipQq- The master equations could be 
derived from (|25p with techniques from non-equilibrium 
physics [36) . In lowest order perturbation theory, the tran- 
sition rates appearing in the master equation would be 
given by (f38|) . In the following, we will not use this ad- 
vanced approach. Instead we will calculate s e and r e per- 
turbatively by appropriately summing and weighting the 
transition rate (1381) over initial and final states. 



3.3 Sticking coefficient 

In order to calculate the sticking coefficient s e we consider 
the positive half space (z > 0) as a kind of quantum- 
mechanical boundary layer (see Fig. [5]) . A measure of the 
tendency SQ n ^Q/ q i with which an electron approaching in 
an unbound state ^Q'q' the plasma-boundary at z = 
gets stuck in a bound state is then the time it takes 
the electron to traverse the boundary layer forwards and 
backwards divided by the time it takes the electron to 
make a transition from ^Q' q ' to <?Q n [38] . 

Since the width of the quantum-mechanical boundary 
layer is L, 



S, 



2L 



1 



Qn.Q' q' 



(QYI^rlQV) 



W- 1 (Qn,Q'q') 



,(46) 



where the denominator in the first factor is the velocity 
matrix element calculated with the incoming part of the 
state ^Q'qi ; n is the normal vector of the boundary point- 
ing towards the plasma and p = —ih\7, the quantum- 
mechanical momentum operator. Using the asymptotic 
form of the unbounded wavefunctions given in Eq. (|85[) 
of appendix B, wc find 



Hence, 



Sc 



i(QVI— IQVh 



16Lm e aB 



Hq' 



hq' 



8m e a,B 



W(Qn,Q'q') 



(47) 



(48) 



The tendency with which the electron approaching the 
boundary in the state gets stuck in any one of the 

bound states - the energy resolved sticking coefficient - is 
then simply given by 



SQ'q' 



^SQn^'q' 
Qn 

16Lm e aB 



hq' 



J2w(Qn,Q'q') . 



(49) 



The sticking coefficient s e entering the rate equation ^ 
is an energy-averaged sticking coefficient resulting from 
an appropriately performed sum over Sq'q' . As mentioned 
before, a rigorous derivation of an expression for s e should 
be based on the master equation for the occupancies of the 
surface states [351136] . A simpler way to obtain s e is how- 
ever to regard the wall as a particle detector. The global 
sticking coefficient can then be defined as 



Q'q> 



SQ'q'q'riQ 



Q'q' 



(50) 



where riQi q i are the occupancies of the unbound surface 
states tyQ'qi ■ 

The occupancies riQi q i depend on the properties of the 
plasma. It is tempting to simply identify riQi q i with the 
incoming part of the electron distribution function as it 
arises on the surface from the solution of the Boltzmann- 
Poisson equations. However, one should keep in mind that 
the distribution function is a classical object whereas riQi q i 
is a quantum-mechanical expectation value. There arises 
therefore the question how the quantum-mechanical pro- 
cesses encoded in the above equations can be properly fed 
into the semiclassical description of the plasma in terms of 
Boltzmann-Poisson equations. The issue is subtle because 
at the plasma boundary the potential varies so rapidly 
that the basic assumptions of the validity of the Boltz- 
mann equation no longer hold. Mathematically, the mi- 
crophysics should be put into a surface scattering kernel, 
course-grained over a few nm, which connects, generally 
retarded in time, the incoming electron distribution func- 
tion with the outgoing one. But even for neutral particles, 
a microscopic derivation of such a scattering kernel has 
not yet been given. There exist only more or less plausi- 
ble phenomenological expressions which parameterize the 
kernel with accommodation coefficients [75 . 

From the boundary-layer point of view used in the 
derivation of Eqs. (|46|) -([50 |l . the plasma, or, more pre- 
cisely, the sheath of the plasma, is infinitely far away from 
the plasma boundary. Rigorously speaking, we can thus 
say nothing about how the microphysics at the plasma 
boundary merges with the physics in the plasma sheath. 

To make nevertheless contact with the plasma we have 
to guess how the unbound surface states ^Q' q > are oc- 
cupied. For simplicity we assume Maxwellian occupancy, 
with an electron temperature T e = (fcs/3 e ) _1 , but other 
guesses, more appropriate for the plasma sheath, are also 
conceivable. For Maxwellian electrons, the global sticking 
coefficient is given by 



Eqv S Qlq ,q' exp[-/3 e £ QV ] 
Eqv Q' exp[-/3 e £ QV ] 



(51) 



nQ 



In the limit L —> oo and A — > oo the momentum sum- 
mations in Eqs. (146")) - ([ST]) become integrals. The calcula- 
tion of SQi q i and s e reduces therefore to the calculation 
of high-dimensional integrals. In appendix C we describe 
the approximations invoked for the integrals. Some of the 
integrals can then be analytically performed. But the fi- 
nal expressions for the sticking coefficients remain multi- 
dimensional integrals, which have to be done numerically. 
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Energy [eV] 

Fig. 9. Energy resolved sticking coefficient for an electron hit- 
ting perpendicularly an aluminum surface at fcsT s = 0.05eV. 
The screening wavenumber for the Coulomb interaction be- 
tween an incident plasma electron and an internal aluminum 
electron, (fc s ) sur f acc , is not well known. Results are therefore 
shown for (fc s ) sur f ace /(fc s )buik = 0.1 (weak screening, strong 
coupling) and for (fc s ) aur f aC e/(fc s )buik = 0.6 (moderate screen- 
ing, weak coupling); (fc s )buik is the screening wavenumber of 
aluminum (see table [2}. Since (fc s ) surfaco = 0.6(fc a ) bu ik is most 
probably the relevant screening parameter 37,40 , the sticking 
coefficient is rather small. 

Measuring energies in units of Ro and distances in 
units of as, Eq. (|5T|) for the global sticking coefficient 
reduces to 

/4\ 2 /3 3/2 
s * ={-) ^ W , (52) 

where 

rOO r OC 1 I ( \ 

/ dR du 2 ,*,[ l 2 h{R,u)g(R^) (53) 

JO J -co 1 + l-fty K s) 

with tib(E) — l/(exp[/3 s _E] — 1) the Bose distribution 
function and h(R,u>) and g(R,u>) two functions defined, 
respectively, in appendix C by Eq. (|104[) and (|105|) . 

Below we also present results for the energy resolved 
sticking coefficient for perpendicular incidence (Q' — 0). 
It is given by 

^' = (-) ^ J Q dRg (R,E') , (54) 

with E' — q' 2 and g ± (R,E') a function defined in ap- 
pendix C, Eq. (fTU7|) . 

The functions h(R,uj), g(R,uj), and g ± (R,E') con- 
tain summations over the Rydberg series of bound surface 
states. If not stated otherwise, we truncated these sums 
after N — 15 terms. These functions are moreover de- 
fined in terms of integrals which can be done only numeri- 
cally. We use Gaussian integration with 40—80 integration 
points. More specifically, h(R,tu) and g ± (R,E') are one- 
dimensional integrals while g(R,u) is a two-dimensional 
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Fig. 10. Dependence of the energy resolved sticking coefficient 
for perpendicular incidence on the number N of bound states 
included in the calculation. Except of the screening wavenum- 
ber, which is set to (fc s ) sur f aco = 0.6(fc s )buik, the parameters are 
identical to the ones used in Fig. [9] The dashed lines are for 
N — 10 and N = 20, respectively, indicating the fast converge 
with respect to N. 



one. Hence, s e and Sjjj, are given by a five-dimensional 
and a two-dimensional integral, respectively. 

In the formulae for the sticking coefficients we mul- 
tiplied the binding energies of the surface states \E n \ ob- 
tained from Eq. (f2"5|) by an overall factor of 0.7. This value 
was chosen to bring the binding energy of the lowest sur- 
face state \Ei\ = 0.85ey in accordance with the experi- 
mentally measured value for copper: |i?i| Cu 0.6eV [5"T] . 
For other metals we used the same correction factor. 

Figure [9] shows the results for S^, when an electron 
with energy E' hits perpendicularly an aluminum bound- 
ary at fcflTs = 0.05ey. Representative for weak and mod- 
erate screening we plotted data for (fc s ) sur f a ce/(&s)buik 
1). 1 and (A: s ) sur f aco /(fc s )buik = 0.6. The latter is the screen- 
ing parameter used in [BTlPfO"] to study the interaction of 
positrons with an aluminum surface. If the corresponding 
value for l/(fc s ) SU rface is indeed a reasonable estimate for 
the length on which the Coulomb interaction between an 
external and an internal electron is screened, the sticking 
coefficient for electrons should be extremely small, of the 
order of 10 -4 . Only for weak screening, and thus strong 
coupling, does Sjjj, approach values of the order of 10 _1 
which are perhaps closer to the value one would expect on 
first sight. 

To clarify the contribution the various bound states 
have to the sticking coefficient, we plot in Fig. [10] the de- 
pendence of S^, on the number N of bound states included 
in the calculation. As can be seen, the lowest bound state 
(N = 1) contributes only roughly 40% to the total S^,. 
The sticking coefficient increases then with increasing N 
but converges for N as 10 — 20. Because of this fast con- 
vergence we present all results below only for N — 15. 
The reason for the convergence can be traced back to the 
decrease of the electronic matrix element (Q) defined 
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Fig. 11. The global sticking coefficient s e for a thermal beam 
of electrons with ksT^ = 5eV hitting various metal surfaces 
at UbTs = 0.05eV as a function of (fc s ) sur f ace /(fe s )b u ik, where 
(fcs)buik is the screening wavenumber in the bulk of the respec- 
tive metal (see table [2}. Following |37H40j . we would expect 
(fcs)surface = 0.6(fc s )b u ik to be a reasonable estimate for the 
screening parameter. Hence, s e ~ 10 -5 — 10~ 4 . 



eigenenergies and eigenstates of the 1/z potential and not 
the ones of an artificial box potential. 

Usually it is assumed that s e is also at least of the 
order of 0.1 [65] . This expectation seems to be primarily 
based on the semiclassical back-on-the envelop-estimate 
of Umebayashi and Nakano [76] . It is thus appropriate to 
discuss their approach in some detail. 

From the energy AE S an electron can exchange in a 
single classical collision with the constituents of the solid 
they first estimated, using the analogy to the Mossbauer 
effect, the probability a for inelastic one-phonon emis- 
sion. For that purpose, they had to estimate the num- 
ber N c of constituents of the surface an electron with 
a de-Broglie wavelength corresponding to its kinetic en- 
ergy Eq, Ag B = 2TraB\/Ro/E , simultaneously impacts. 
A rough estimate is N c = (Af s /a) 2 , where a is the lattice 
constant of the material. Under the assumption that the 
electron hops along the surface they then calculated the 
probability with which the electron does not escape after 
I hops where I is the number of inelastic collisions which 
are necessary for the electron to transfer its whole positive 
kinetic energy to the lattice, that is, to end up in a state 
of negative energy. Identifying this probability with the 
(global) sticking coefficient, they obtained 



in Eq. (|44| . which we approximate by I^x n {Q = 0) (see 
appendix C), with increasing n, where n — 1,2, ... labels 
the bound surface states. 

Global sticking coefficients s e as a function of the screen- 
ing wavenumber (fc s ) sur f aC c are shown in Fig. QT] for dif- 
ferent metals. For (fc s ) sur f ace /(fc s )buik > 0.4, the sticking 
coefficients are again extremely small. As expected they 
increase with decreasing (fc s ) sur f acc /(fc s )buik, reaching val- 
ues close to unity for weak screening. In this strong cou- 
pling regime, our perturbative calculation of s e is no longer 
valid. We believe however that (fc s ) B urface/(fc s )bulk < 0.4 
is unphysical. The kink around (fe s ) sur face/(^s)bulk ~ 0.25 
must be due to an accidental resonance in g(R, cu). It is of 
no physical significance. 

Why is the sticking coefficient for electrons so small? 
We have no satisfying explanation. Our calculation pro- 
duces small a sticking coefficient because the matrix ele- 
ment Q37[) turns out to be very small. We certainly under- 
estimate it because the wavefunctions of the approximate 
potentials (|39|) and HI]) vanish in complementary halfs- 
paces, in contrast to the exact wavefunctions which have 
tails. Nevertheless it is hard to image the tails of the wave- 
functions to increase the matrix elements by three orders 
of magnitude. 

The approximations we had to make to end up with 
manageable equations for s e , in particular, the assump- 
tions about the momentum dependence of the electronic 
matrix elements (see appendix C and, for a discussion, the 
next section) should also not lead to a sticking coefficient 
which is more than one order of magnitude off. In this 
respect let us emphasize that in contrast to the calcula- 
tions performed in for a positron, which produce 
positron sticking coefficients of the order of 0.1, we use the 



i-i 

where $ = (E — iAE)/Eb is the escape probability af- 
ter i inelastic collisions [77], AE = 2AE s /3N c a, AE S = 
4m e (Eo + Eb)/M, E^ is the depth of the surface potential 
and M is the mass of the constituents of the solid. 

Sticking coefficients for graphite obtained from Eq. (|55|) 
are shown in Fig. [T2] Within Umebayashi and Nakano's 
semiclassical approach we identified Eb with the binding 
energy of the electron. According to Fig. [T^] the sticking 
coefficient very quickly approaches extremely small val- 
ues with increasing energy Eq. The smaller the binding 
energy Ef,, the faster the decrease. The values for s e orig- 
inally given by Umebayashi and Nakano were for kinetic 
energies smaller than 0.0026eU and binding energies larger 
than leV. Only in this parameter regime is the sticking 
coefficient close to one. In the parameter range which is 
of interest to us (kinetic and binding energies at least a 
few tenth of an electron volt) Umebayashi and Nakano's 
estimate gives also an extremely small sticking coefficient. 

We should of course not directly compare the results 
obtained from Eq. ([52]) with the ones obtained from Eq. (|55|) 
because Eq. (|52[) assumes energy relaxation due to internal 
electron-hole pairs whereas Eq. (|55|) assumes energy re- 
laxation due to phonons. However, a quantum-mechanical 
calculation of the phonon-induced electron sticking coeffi- 
cient at vanishing lattice temperature also shows that s e w 
10 -4 [52], in contrast to what Umebayashi and Nakano 
find. Although they incorporate some quantum mechanics 
their approach is basically classical. It is based on the no- 
tion of a classical particle hopping around on the surface 
and exchanging energy with the solid in binary encoun- 
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Fig. 12. Electron sticking coefficient obtained from Ume- 
bayashi and Nakano's phenomenological model [76], see 
Eq. (|55p . The solid lines are for the e~:graphite system origi- 
nally considered by them (ksTo = 420K, Mc = 12m p , where 
m p is the proton mass, and a = 2.5A) and the dashed line is for 
an e":Cu system (E b = 0.6eV, k B T D = 343^, M Cu = 64m p , 
and a = 3. 61 A). The sticking coefficient diminishes rapidly 
with increasing electron energy and approaches one at zero 
electron energy, in contrast to what one would expect from a 
quantum- mechanical calculation [39] . 



ters. As in any classical theory for the sticking coefficient, 
it is therefore not surprising that the sticking coefficient 
they obtain approaches unity for the low energies they 
consider [55] , 



3.4 Desorption time 

We now calculate the electron desorption time r e . For that 
purpose, we have to specify the occupancies of the bound 
electron surface states. In general, this is a critical issue. 
However, provided the desorption time r e is much larger 
than the time it takes to establish thermal equilibrium 
with the boundary, it is plausible to assume that bound 
electron surface states are populated according to 



n Qn ~ e-Kp[-l3 s E Qn } 



(56) 



where T 8 — l/kgf3 s is the surface temperature. 

Desorption is accomplished as soon as the electron is in 
any one of the unbound surface states. Hence, the inverse 
of the desorption time, that is, the desorption rate, is given 
by [36] 



1 = EQ'n'E Qq ^P[-PsE Qln/ }W(Qq,Q'n') 



(57) 



where W(Qq, Q'n') is the transition rate from the bound 
surface state (Q',n') to the unbound surface state (Q,q) 
as defined by the golden rule ([3"5[) . 

Measuring again energies in units of Rq and distances 
in units of as and using the same approximations as in 
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Fig. 13. Desorption time r e for an electron bound in the 
polarization-induced external surface states of various metal 
surfaces at ksT s = 0.05eF as a function of (fc s ) sur f aC e/(fcs)buik- 

Since (&s)surfacc = 0.6(/c s 

)buik is most probably the relevant 
screening wavenumber [371140) , r e ~ 10 _2 s. For the used mate- 
rial parameters, see table [2] 



the calculation of the sticking coefficient (see appendix C) 
the desorption rate can be cast into 



2tt 3 hZ 



dcsorb 



(58) 



where 



/dcsorb = / dRl du * * f(R, u)g(R, u) ,(59) 

Z = J2 n exp[— /3 s E n ], ub{E) is again the Bose distribu- 
tion function, and f{R,u>) is an one-dimensional integral 
defined in appendix C, Eq. (jllip . Thus, to obtain r^ 1 from 



Eq. ([58]) we have to do a five-dimensional integral. As for 
the calculation of s e we again use Gaussian quadratures 
for that purpose. 

In Figure [131 we present, as a function of the screening 
parameter, numerical results for r e for an electron bound 
in the polarization-induced external surface states of var- 
ious metal surfaces at fcflT s = 0.056^. To be close to 
reality, we again corrected the binding energies \E n \ by 
a factor 0.7. As can be seen, except for small screening 
parameters and thus strong coupling, r e w lCP 2 s. 

Compared to typical desorption times for neutral mole- 
cules, which are of the order of 10 _6 s or less [36], the 
electron desorption time we find is rather long. This is a 
consequence of our assumption that the bound electron 
is in thermal equilibrium with the surface (viz: Eq. (j56[l ) 
and the fact that the binding energy of the lowest surface 
state | .Eil 3> kgTg. Thus, the electron desorbs de facto 
from the lowest surface state which has a binding energy 
of ~ O^eT^. The binding energies for neutral molecules, on 
the other hand, are typically one order smaller and thus 
of the order of ksT s resulting in much larger desorption 
rates and thus shorter desorption times. 
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Fig. 14. The product s e r e as a function of fcsT e for a ther- 
mal beam of electrons hitting a copper surface at various 
temperatures fc_gT s . The surface screening wavenumber is set 
to (fc s ) surfacc = 0.6(fc s ) bu i k , where (fc s ) bu ik is the screening 
wavenumber for copper (see table [2}. 



In the model for the quasi-stationary charge of a dust 
particle presented in the previous section the product (sr) e 
was of central importance. Combining (f52|) and (|58|) . the 
microscopic approach gives 



/ \ 11 -LU-'stick 

[ST)S = (fc B T e )3/2(fc B T s )-l/2 I dcaolh 



exp[/3 s |£: n |](60) 



where /stick and Idesorb are defined in Eqs. |53|) and (|59|) . 
respectively. 

Figure fT4l shows numerical results for (sr) e for a copper 
surface as a function of the electron and surface tempera- 
ture. The screening wavenumber (fc s ) SU rface = 0.6(fc s )buik 
and the binding energies are again corrected by the fac- 
tor 0.7 which makes \E\\ to coincide with the experimen- 
tal value for copper. Notice, the weak dependence of the 
product (sr) e on the electron temperature and the rather 
strong dependence on the surface temperature. The latter 
is of course a consequence of the exponential function in 
Eq. ([60| . Although the sticking coefficient and desorption 
times have values which are perhaps in contradiction to 
naive expectations, s e being extremely small and r e being 
rather large, the product (sr) e has the order of magnitude 
expected from our surface model (see section 2). In par- 
ticular, (sr) e ~ 10~ 6 s for A:bT s = 0.045ey would produce 
grain charges of the correct order of magnitude. Thus, us- 
ing Eq. ([60)1 instead of Eq. (fl3)) and ksT s as an adjustable 
parameter, which is still necessary because the grain tem- 
perature is unknown, we could produce, for physically re- 
alistic surface temperatures, surface charges for metallic 
grains which are in accordance with experiment [19j . 

Although the microscopic Eq. (f6TJ)) has a similar struc- 
ture as the phenomenological expression (|13|) there are 
significant differences. First, the microscopic formula con- 
tains more than one bound state and depends not only 
on T, but also on T e . In addition, there is a numerical 



factor 16 = 2x8 where the factor 2 comes from the 
fact that an electron traversing the quantum-mechanical 
boundary layer can make a transition to a bound state on 
its way towards the surface and on its way back to the 
plasma and the factor 8 originates from the asymptotic 
form of the wavefunction for the incoming electron. The 
phenomenological approach simply assumes here a plane 
wave whereas the microscopic approach works with Whit- 
taker functions (see appendix B and C). Most importantly, 
however, the two functions J s tick and idesorb, which depend 
on the microscopic details of the inelastic scattering pro- 
cess driving physisorption, and thus on the electron and 
surface temperature as well as material parameters such 
as the screening wavenumber, are in general not identical. 

Hence, /stick /^desorb 7^ 1- 

For the hypothetical case of a single bound state, how- 
ever, whose binding energy \E\\ is much larger than ksT s 
and ksT e , Eq. (frjO"]) reduces to a form which, for ksT,, = 
ksT s , becomes identical to the phenomenological expres- 
sion (|13p . except of the numerical factor referred to in the 
previous paragraph. The simplification arises because for 
low temperatures the integrals defining / s tick and Idesorb 
can be calculated asymptotically within Laplace's approx- 
imation (see appendix C). The sticking coefficient and des- 
orption time are then given by 



7T/V Pe 

L a 2„2 h eXp[0 s \Ei\] 
T e = 87T fj s 



i?,n 



r(!)| 2 ; 



(61) 
(62) 



with g defined in appendix C, Eq. (|1 18[) . Thus, the prod- 
uct, 



(sr) 



16h 



(k B T s )^(k B T e )-^ 



exp^rl] , (63) 



is independent of the microscopic details of the inelastic 

scattering processes encoded in the product |J^| 2 5. Iden- 
tifying \Ei\ with the electron desorption energy Ef and 
setting ksT e — ksT s , we finally obtain, except of the nu- 
merical factor 16 = 8 x 2, from Eq. ([63]) the phenomeno- 
logical expression (Tl"3"]) . 

Using Eqs. (f6"Tj) - (|63[) we find for an electron at a copper 
boundary with k B T e = k B T s = 0.045eU, s L e = 6.23xl0~ 6 , 
t c l = 0.131s, and (sr)£ = 8.17 x 10~ 7 s. Taking only one 
bound state into account, the corresponding values ob- 
tained from Eqs. $52$ and ([55]) are s e = 4.42 x 10~ 6 and 
r e = 0.135s, which leads to (sr) e = 6 x 10~ 7 s, indicating 
that at low temperatures Laplace's approximation works 
indeed reasonably well. Since r e does not depend on ksTe 
and ksT s is usually much smaller than \E\\, approxima- 
tion (|62| for r e can be actually always applied, provided 
the assumption is correct, that the electron is initially in 
thermal equilibrium with the surface and hence basically 
in its lowest bound state. The approximation (|6ip for s e , 
on the other hand, deteriorates quickly with increasing 
electron temperature, as does the approximation (|63[) for 
(sr) e . 



18 



F. X. Bronold et al.: Physisorption kinetics of electrons at plasma boundaries 



It is reassuring to be able to derive, under certain con- 
ditions and except of a numerical factor, whose origin is 
however clear, from the microscopic expressions for s e and 
r e the phenomenological relation (fT"3| we used in [29] as 
an estimate for (sr) e - That ([50)1 can be reduced to (fT3"]l 
is a consequence of the perturbative calculation of s e and 
r e using the golden rule transition rate (|38|) which obeys 
detailed balance. In this respect, our calculation is on par 
with Lcnnard- Jones and Devonshire's original microscopic 
derivation of the product sr for a neutral adsorbant [30] . 
In contrast to them, we keep however the lateral motion 
of the adsorbing particle and the elementary excitations 
of the solid responsible for energy relaxation are electron- 
hole pairs and not phonons. 



3.5 Critique 

In the previous subsections we demonstrated for a partic- 
ular CclS6, cl metallic boundary with electron energy relax- 
ation due to creation and annihilation of internal electron- 
hole pairs, how a quantum-mechanical calculation can be 
set up to obtain s e and r e from a microscopic model for the 
electron- wall interaction. To obtain manageable equations 
we had to make various approximations, some were purely 
technical, but others concerned the physics. We now re- 
state and criticize the approximations in the hope that it 
will be read as a list of to-do's. 

We start with the purely technical approximations. In 
the calculation of the transition rate we neglected the de- 
pendence of the matrix element (|37p on the lateral momen- 
tum transfer and approximated furthermore ij^J (Q = 0) 
by its leading term for fc<l. Both approximations can be 
avoided but the final equations become more complex and 
the costs for their numerical handling accordingly higher. 
At the present stage of the investigation this seemed to 
us not justified. Even more so because we do not believe 
that these approximations are the cause for the unexpect- 
edly small values for s e and the unexpectedly large val- 
ues for r e . Neglecting the dependence on the lateral mo- 
mentum overestimates even the matrix elements, hence 
the transition rate, and thus, eventually, s e and t" 1 . The 

fc— dependence of ity (see Eq. ([87]) in appendix B), on the 
other hand, can also not be so large that it increases the 
transition rate by three orders of magnitude as it would 
be required to obtain s e ~ 0.1 — 1 and r e ~ 10~ 5 — 10 _6 s, 
the values one would perhaps naively expect. 

More critical for the matrix element (f37|) are the re- 
placements (|39|) and ([41"]) because they lead to wavefunc- 
tions vanishing in complementary halfspaces and thus to 
the factorization ([43]) of the matrix element. In reality the 
wavefunctions have tails in the complementary halfspaces. 
A model neglecting the tails underestimates therefore the 
matrix element. In addition, the replacements lead to the 
loss of crystal-induced surface states, about which we have 
more to say below, and hard-wire the artificial treatment 
of surface and volume electrons as two separate species. A 
more realistic modeling should therefore avoid these two 
approximations . 



Both the electron sticking coefficient s e and the elec- 
tron desorption time r e were obtained from the golden 
rule for transitions between bound and unbound surface 
states. This is only justified for weak coupling and when 
one quanta of elementary excitation suffices for the transi- 
tion. When the coupling is strong, or when more than one 
quanta are necessary, a generalized golden rule has to be 
used in which the interaction matrix element is replaced 
by the corresponding on-shell T- matrix 3 1 , 36J . The calcu- 
lation becomes more tedious but it can be done. A prin- 
ciple shortcoming, however, of any approach which uses 
golden-rule-type transition rates directly to calculate s e 
and r e is that it assumes the occupancies of surface states 
to be only weakly affected by the transitions itself. From 
the physisorption of neutral particles it is known that this 
is in general not true [34] . 

The calculation of the desorption time, for instance, 
was based on the assumption that the desorbing elec- 
tron is initially in thermal equilibrium with the surface 
and that during the desorption process the equilibrium 
occupancy of the surface states does not change. The des- 
orption time is thus much larger than the timescale on 
which thermal equilibrium at the surface is established, 
in which case the electron basically always desorbs from 
the lowest bound surface state. The equilibration on the 
surface is controlled by transitions between bound surface 
states. They have to be much faster than transitions be- 
tween bound and unbound surface states. In the golden 
rule approach this information is put in by hand. Thus, 
although the r e obtained is consistent with the equilibrium 
assumption, it does not justify it. For that purpose, the 
calculation of r e has to be based on quantum-kinetic mas- 
ter equations which include not only transitions between 
bound and unbound surface states but also transitions be- 
tween two bound surface states [3~5|J36j . 

Master equations are also required when the elemen- 
tary excitations of the solid have not enough energy to 
couple the lowest bound surface states to the continuum. 
In that case, the cascade model developed by Gortel and 
coworkers [35] has to be used. Its main idea is that an elec- 
tron initially bound in a deep state can successively climb 
up to the continuum using weaker bound states as inter- 
mediaries. By necessity, it thus also contains transitions 
between bound surface states. 

For metals internal electron-hole pairs provide the most 
efficient electron energy relaxation channel, with phonons 
and other elementary excitations being unimportant, be- 
cause their energy is either too high (plasmons) or too low 
(phonons). Both leads to severe restrictions in the avail- 
able phase space. For dielectric boundaries, however, it is 
the energy of internal electron-hole pairs, whose energy is 
of the order of the intrinsic energy gap, which is too high 
for having any effect. Electron energy relaxation should 
then be primarily driven by phonons. Their energy, how- 
ever, is in the cases of interest, for instance, graphite or 
silicon, too low for promoting an electron from the low- 
est surface states all the way up to the continuum. Hence, 
s e and r e have to be calculated from Gortel et aVs cas- 
cade model. When the energy of the phonon is moreover 
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not enough to connect two neighboring bound states, the 
transition rates entering the master equation have to be 
obtained from the generalized golden rule containing the 
T-matrix for electron-phonon coupling. 

We expect multiphonon processes to play a role for all 
dielectric boundaries even for graphite boundaries, where 
the Debye energy is rather high, but not high enough to 
couple the two lowest image states (see table [J). This cou- 
pling, on the other hand, is the rate-limiting one, that is, 
the one which determines the electron desorption time. 
Multiphonon processes remain important when in addi- 
tion to image states also crystal-induced surface states or 
dangling bonds are included because these states, being 
stronger bound than image states, are energetically deep 
in the gap and thus far away from the vacuum (plasma) 
level. 

In our model we made the overall assumption that 
plasma electrons cannot enter the plasma boundary (hard 
boundary condition at z = 0). At least electrons with an 
energy larger than the projected energy gap of the solid, 
can however enter the plasma boundary, scatter inside the 
material, before bouncing back to the surface, where they 
may be either re-emitted to the plasma or trapped in sur- 
face bound states. Processes inside the material can be 
thus only neglected when the projected energy gap is much 
larger than the typical energies of plasma electrons, that 
is, Eg ^> fcsT e , and when the floating potential U is ap- 
proximately in the middle of this large gap. 

Here we come to a potentially very interesting point, 
in particular, as far as metallic plasma boundaries are 
concerned. According to Fig. [5] the projected energy gap 
depends on the crystallographic orientation of the sur- 
face. Even planar metallic plasma boundaries will how- 
ever almost never coincide with a single crystallographic 
plane. At best, they contain large, crystallographically 
well-defined facets, as we discussed in the context of spher- 
ical grains. Hence, the projected energy gap varies along 
the boundary. Regions can thus be expected where surface 
states are absent and plasma electrons can easily enter the 
boundary. In other regions large gaps prevent plasma elec- 
trons from entering the material. Instead they would sit in 
surface states. How all this affects the spatial distribution 
of surface charges is an open question. 

For dielectric boundaries the projected energy gap is 
of the order of the intrinsic gap of the bulk material. It 
depends only weakly on the crystallographic plane. But a 
problem which concerns both metallic and dielectric sur- 
faces is the existence of crystal-induced surface states. We 
would expect them to be less important for physisorp- 
tion of electrons. Being strongly bound and having a cen- 
ter of gravity very close to the surface or even inside the 
plasma boundary, the spatial overlap between unbound 
surface states and crystal-induced surface states should 
be rather small. Hence, the matrix element controlling 
sticking into or desorption from a crystal-induced surface 
state should be much smaller than the corresponding ma- 
trix element involving weakly-bound polarization-induced 
surface states, which are always exterior to the boundary 
and, on a microscopic scale, even relatively far away from 



the surface. However, only a detailed study can show if 
our intuition is correct. 

Another problem concerning both metallic and dielec- 
tric plasma boundaries is surface roughness. In our model, 
the plasma boundary is a well-defined mathematical plane. 
On the atomistic scale, and we actually do calculations on 
this scale, the surface is however not perfect. In a refined 
model for surface states this aspect, possibly in conjunc- 
tion with surface reconstructions El and chemical contam- 
ination has to be taken into account. 

Throughout we implicitly assumed that bound sur- 
face states exist although the exact surface potential sup- 
porting them is unknown. Since surface states have been 
detected many times [3|45llli|4^^ 
[55] this assumption seems to be justified. Naturally, it 
would be desirable to calculate the surface potential from 
first principles. However, if not illusionary, it is at least 
very challenging, even when the plasma boundary is pla- 
nar and crystallographically well defined. The quantum- 
mechanical exchange and correlation effects determining 
the tail of the surface potential are beyond the local- 
density approximation, the work-horse of most ab-initio 
packages for the calculation of the electronic structure of 
solids. Instead, non-local density functional theory [73] has 
to be used which is much more complicated to implement. 
A compromise would be to calculate the potential inside 
the boundary from an ab-initio local-density package and 
then continuously match this "internal" potential to the 
"external" potential deduced from a model Hamiltonian 
of the type presented in appendix A. The model produces 
a diverging potential only in the simplest approximation. 
With methods adapted from bulk polaron theory 42.43 
potentials could be deduced which are finite at the bound- 
ary and thus continuously matchable with the periodic 
crystal potential obtained from the local density approxi- 
mation. 

As in other branches of surface science [361167]. a gen- 
eral strategy to short-circuit unknown microscopic details 
about the surface would be to work with simple, possibly 
analytically solvable models containing parameters that 
can be adjusted to experimentally measured quantities, 
for instance, the binding energy of surface states. 

For this strategy to work, experimental techniques suit- 
able for directly probing the electronic properties of sur- 
faces, for instance, inverse photoemission spectroscopy [44, 
45,47], from which the binding energy and the lifetime 
of unoccupied electron surface states can be determined, 
have to be adapted to plasma boundaries. In addition, 
macroscopic quantities, such as the quasi-stationary sur- 
face charge, the surface temperature, and the tempera- 
ture and density of the electrons in the plasma have to 
be also known. So far, however, these combined data are 
not available for any experiment. For sure, surface charges 
have been measured in dielectric barrier discharges :25, 24, 
23,22 and of course in complex plasmas, where in fact a 

4 Here we do not mean the reconstruction of the surface due 
to impacting plasma particles but the intrinsic reconstruction 
leading to geometrical differences between real terminations of 
crystals and ideal crystallographic planes [67] . 
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great variety of techniques has been invoked to determine 
the charge of floating /im-sized dust particles p"3jlT4, 15. 
[T6lfT7lfT8l[T9] . But particularly in the experiments measur- 
ing grain charges the diagnostics of the hosting plasma 
is usually missing. In addition, although it is possible to 
measure the temperature of the grain [72] , grain temper- 
ature and charge have not yet been measured simultane- 
ously. For the microscopic modeling of surface charges it is 
however important to know at least these two quantities. 



4 Concluding remarks 

In this colloquium we proposed to treat the interaction 
of electrons and ions with inert plasma boundaries, that 
is, boundaries which stay intact during their exposure to 
the plasma, as a physisorption process involving surface 
states. The sticking coefficients s e ,i and desorption times 
T et i can then be calculated from microscopic models con- 
taining (i) a static potential supporting bound and un- 
bound surface states and (ii) a coupling of these states to 
an environment which triggers transitions between them. 
Microscopically, the sticking of an electron or ion to the 
surface corresponds then to a transition from an unbound 
surface state to a bound one. Desorption of an electron or 
ion from the wall is then simply the reverse process. 

Although this point of view can be applied to ions 
and electrons, we worked it out - for the particular case 
of a metallic boundary and within the simplest possible 
model - only for electrons because the surface states for 
electrons are surface states in the ordinary sense, that is, 
states which are only a few nanometers away from the 
surface. The environment responsible for transitions be- 
tween electron surface states, and thus for sticking and 
desorption of an electron, are therefore the elementary 
excitations of the solid. For ions, however, as soon as the 
surface collected some electrons, the surface potential is 
the long-range attractive Coulomb potential. Sticking and 
desorption of ions occurs thus far away from the surface. 
Nevertheless, provided the surrounding plasma is taken as 
the environment triggering transitions between ion surface 
states, the dynamics and kinetics of ions in front of the 
boundary can be described in close analogy to the elec- 
tron dynamics and kinetics occurring much closer to the 
surface. Since without the surface no attractive Coulomb 
potential for ions would exist, the ion dynamics and ki- 
netics is also a kind of surface physics although it takes 
place far away from the surface. 

Ions are much heavier than electrons and the potential 
most relevant for them, the Coulomb potential, varies on 
a scale much larger than the ion's de-Broglie wavelength. 
Quantum mechanics is thus not really required for study- 
ing the ion kinetics in front of a plasma boundary. Instead 
of pushing the quantum-mechanical techniques we used for 
electrons to the semiclassical regime, it is thus also pos- 
sible to analyze ions with Boltzmann equations. In that 
case it is however crucial to set up two Boltzmann equa- 
tions, one for unbound ions and one for bound ions. As in 
the quantum-mechanical calculation, collisions of bound 



and unbound ions with the atoms/molecules of the back- 
ground gas determine the number of trapped ions and how 
they are spatially distributed. 

Studying the ion dynamics and kinetics is important 
because it affects the rate with which ions and electrons 
may recombine in the vicinity of the grain surface. If the 
corresponding flux an<j e (Ji is larger than the electron des- 
orption flux t~ a e , the charge of the grain is the one which 
balances on the grain surface the electron collection flux 
s e jP lasma with the ion collection flux Sjjf lasma and not 
with the electron desorption flux t~ a e as in our surface 
model for the grain charge. Provided s e ~ s$ this would 
eventually lead to the standard criterion from which the 
grain charge is calculated. We emphasize in this respect 
however that the rate equations for the surface densi- 
ties a Ey i are phcnomcnological. They should be derived 
from a quantum-mechanical surface scattering kernel tak- 
ing bound surface states into account. Only then would 
we know if the microscopically obtained s e and r e and the 
macroscopic plasma fluxes j^^ sma are as simply connected 
as in the phenomenological rate equations. In any case, 
the quantum-mechanical approach for calculating s e and 
r e stands by itself irrespective of the fate of our surface 
model for the grain charge. 

Admittingly, the microphysics at the plasma-boundary 
we discuss is not the one utilized in plasma technology. 
Precisely the processes we excluded are most important 
there: Implantation of heavy particles, reconstruction or 
destruction of the surface due to high energy particles, 
and chemical modification due to radicals, to name just a 
few. The target surfaces are of course charged but, from 
the perspective of plasma technology, the surface charges 
only control the particle fluxes to the surfaces. Properties 
other than their mere existence are of no concern. 

From a microscopic point of view, the technologically 
important surface processes just listed are extremely com- 
plicated. A description of these processes at a level, let say, 
solid state physicists describe superconductivity in bulk 
metals is certainly far from reach. It may even not be re- 
quired for plasma technology to proceed as a business. But 
as in other branches of science, it is the pleasure and duty 
of research driven by curiosity to push the understanding 
of particular processes, technologically relevant or not, to 
an ever increasing level of sophistication. We firmly be- 
lieve, the microphysics at an inert plasma boundary is 
now ready for a truly microscopic understanding. It is our 
hope to have inspired other groups joining us on our jour- 
ney to the microphysics at an inert plasma boundary. In 
particular, however, we hope to have found experimental- 
ists eager to design experiments with well-defined model 
surfaces, which come as close as possible to the idealized 
boundaries theorists have to consider in their calculations, 
and at the same time are accessible to the surface diag- 
nostics used elsewhere in surface science. 

A Microscopic model for the image potential 

In this appendix we discuss a microscopic model which in- 
terprets the image potential in terms of virtual excitation 
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of surface modes [4T | l42 l l43] . The model is applicable to 
metals and dielectrics. 

To be specific, we consider a planar plasma boundary 
in the xy plane putting the plasma in the positive halfs- 
pace defined by z > 0. A convenient starting point for a 
microscopic description of the polarization-induced inter- 
action between an electron and a boundary is the single 
electron Hamiltonian 41,32], 



is the classical image potential arising from virtual exci- 
tation of surface modes and 
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where A is the three-dimensional Laplace operator, a K is 
the creation operator for the polarization-induced surface 
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After the transformation the Hamiltonian reads 
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is a vector potential giving rise to a "minimal-type" 
namic coupling- the last two terms on the rhs of 
between the electron and the surface mode. 

The first two terms on the rhs of (f6"B")) describe an elec- 
tron in a potential. Diagonalizing these two terms, that is, 
using the eigenstates of Eq. (f26|) with V(z) — > V p (z) as a 
basis and ignoring the nonlinear term ~ A 2 we obtain 
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where 



is the coupling function; K = (K x ,K y ) is a two-dimensional 
wavevector, R = (x,y) denotes the projection of the elec- 
tron position onto the surface, whose area is A, and z is 
the distance of the electron from the surface. 

For metals (e = oo), the relevant surface modes are 
surface plasmons with typical energies of a few electron 
volts, for instance, for copper, Huj s w 2eV [5T] , For di- 
electrics (e < oo), the relevant surface modes are opti- 
cal phonons with energies of a few tenth of an electron 
volt, for instance, for graphite, hui s ~ QAZeV when we use 
lo s = LO Ty /(e + l)/2 with e = 12 and uj t = 0A7eV [78] . 

To approximately separate the static from the dynamic 
interaction, we apply to the Hamiltonian (|64|) the unitary 
transformation [43] 
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describes an electron in classical image states. Thus, with- 
out the dynamic coupling to surface modes, H — ► H e , and 
we would have obtained the model we used for the calcu- 
lation of s e and r e . 

Obviously, the dynamic coupling encoded in the last 
term on the rhs of (|71[) renormalizes the classical image 
states. The eigenstates of the full Hamiltonian - the true 
polarization-induced surface states - are not identical to 
the classical image states. The latter should be consid- 
ered as zeroth order (or bare) eigenstates. Better approx- 
imations can be constructed using methods from polaron 
theory 42, 43J . At large enough distances, however, where 
the residual interaction becomes negligibly small, classical 
image states are reasonably good approximations to the 
true polarization-induced surface states. 

Separating the lateral from the vertical motion accord- 
ing to Eqs. (|2"6"]) and (f2"B"]) the matrix element for the dy- 
namic coupling between the electron and the surface mode 
becomes 



G qq ,(Q,K) = 
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In general, G qq i is non-diagonal. It contains intraband 
(<Z = q') an d interband (q ^ q') transitions. The latter 
could in principle affect the physisorption kinetics of elec- 
trons (understood - for the moment - as transitions be- 
tween bound and unbound bare surface states). This hap- 
pens however only when the energy of the surface mode 
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is comparable to fcsT s , where T s is the surface temper- 
ature, as well as comparable to the energy spacing of 
the bare surface states. Transitions between bare surface 
states are then associated with creating or annihilating 
real surface modes in contrast to virtual modes which 
would only renormalize the energies Eq q and the wave- 
functions ipQq. Physisorption would then be triggered by 
other elementary excitations of the solid, for instance, 
phonons and would moreover take place between renor- 
malized surface states. 

Whether the residual interaction directly triggers ph- 
ysisorption of electrons or not depends on the material. 
For dielectric boundaries, for instance, graphite, Huj s ~ 
AE21 (see table [lj. The nondiagonal elements of G qq i 
could thus indeed be important for the physisorption pro- 
cess, especially at high temperatures. For metallic bound- 
aries, however, the energy of the surface plasmon is a few 
electron volts and thus far too high to play any direct 
role in the physisorption process. The dynamic coupling 
to surface plasmons modifies then primarily the proper- 
ties of the surface states in which physisorption takes 
place (energy, wavefunction) . We neglected these modi- 
fications, although at short distances they are not neces- 
sarily small on the scale of the electron binding energy, 
which is the relevant energy scale. We thereby assumed 
that for all distances, not only for large distances, the true 
polarization-induced surface states can be reasonably well 
approximated by classical image states. 
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B Electronic wavefunctions and matrix elements 

In this appendix we summarize the properties of the electronic wavefunctions for the vertical motion of the external 
electron. The results are well-known and the appendix primarily serves the purpose to fix our notation. 

First, we consider bound surface states. Using y = z/2a,Bn, with n = 1,2, ... the quantum number labelling the 
Rydberg series of bound states, the Schrodinger equation (|28| for the vertical motion becomes 



d2 , / x 



7 + - My) = o , 

whose solutions are Whittaker functions [79) . Hence, the wavefunctions which vanish at z — and for z — > oo are 

^„(z) = N n W nA/2 (y) = exp[-y/2]y(-)"- 1 (n - lJIL^fo) , 



(76) 



(77) 

where N n is a normalization constant and i| l 1 2 1 (y) is an associated Laguerre polynomial [79) . The corresponding 
eigenvalues are E n = — Ro/l&n 2 . 

In order to find the normalization constant, we insert the expansion of the Whittaker function [79] . 



where 



W n<1/2 (y) = exp[-y/2}y n J2 a ^ q > 

3=0 



(-)« I>+1)I» 



" g! r(n-q)r(n-q+l) 
with -T(n) the Gamma function, in the normalization integral, 



1 = 



dz\M*t 



Term-by-term integration leads then to 



N„ 



An 3 r(n) 2 aB 



(78) 



(79) 



(80) 



(81) 



which is the defining equation for Af n needed in appendix C. 

For the continuum states, we use y = ikz/2a,B as an independent variable. The Schrodinger equation (|28p can then 
be reduced to (|76|) with n replaced by —ik~ x . The continuum states with energy = Rok 2 /16 which vanish at z = 
are thus given by 79J 



i> k {z) = N k M_ ik -i A/2 (y) 



(82) 



As for any continuum state, to find the normalization constant N k is somewhat tricky. We could normalize ipk(z) 
on the momentum scale but we found it more convenient to use a box-normalization considering the plasma halfspace 
(z > 0) as a slap of width L with L — > oo at the end of the calculation. Thus, N k is determined from the condition 



dz\Mz)f 



(83) 



To do the normalization integral, we utilize the fact that in the limit L — > oo the contribution to the integral 
coming from small z is negligibly small compared to the contribution coming from large z. Hence, we can replace 
in (j83[) tpk(z) by its asymptotic form for large z: 



= N k 



exp[— 7r/2fc] ... exp[— ir/2k + iir] , ., . 

1 y| exp[ifca;/4] + ^ exp[-ikx/4] 



rii + ik- 1 ) 



r(i - ik- 1 ) 



(84) 
(85) 



where we defined in- and outgoing waves which we need in appendix C for the calculation of s e . 
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The normalization constant is then given by 

N k = 



A4 



Lk(l - exp[-27r/fc]) yfh 



(86) 



which also defines A4 needed in appendix C. 

Having appropriately normalized wavefunctions, we can now calculate the electronic matrix element (|44j) . Although 
we could calculate (|44|) for any R and any k we give only the result for R = and k <C 1 because in the calculation 

of s e and r e we eventually approximate (|44|) by -f^L;i(0). The multidimensional integrals defining s e and r e are then 
easier to perform. 

The matrix element we need is 



I^{Q) = 2nM v M k I dyexp[-y/d}W rhl/2 (y)M_ lk -i A/2 (ikny) 



(87) 



Approximating A4 ~ (n/k) 1 / 2 for k <C 1 gives 



7T 1 



with 



h 



dyexp[-y/d]jy„ il/2 (y)M_ lfe -i 4/2 (ifcm/) 



(89) 



which, to be consistent, we also have to calculate for fe<Cl, 

To determine the integral 7i, we use the expansion (j78j) for W n! i/ 2 (y) together with the expansion [79] 

M_ ifc - M/2 (ifcny) = *fcn C m i^) ( ? I +1)/2 y (m+1)/2 J m+ i (2^2/(1 - **0) 

m— L V /J 

for M_ ik -i 1 / 2 (ikny) 1 where C m are constants and J n (y) are Bessel functions 0. Thus, 



(90) 



oc n 



h = ikn } a q C m — — 
z — ' ^—^ nil 

m=0 q=0 1 



(ikn) r 



ife )](m+l)/2 



(91) 



with an integral I 2 which can be found in 



dy exp[-(l/d + l/2)y] exp[n - g + (m + l)/2]J m+1 (2y/ny(l - ik)) 



2d 
2 + d 



n— q-\-m-\-2 



(n-g)![(l-ifc)n] (m+1)/2 



exp 



- (1 - ifc) 



2nd 



2 + d 



L 



(m+l) , 



2nd 



(92) 



with d = l/2nk s . 

Inserting (f9"2")) for k <C 1 into (|9ip and using Cq = 1 [75j we finally obtain 



l4 1 fc ) <i(o)| 2 = fc|4 1) l 2 = ™*( 



2d 
2 + d 



l/nl 



with 



i/«i 2 = E 



(-)« r(n+ 1) / 2d 



9=0 



q! r(n-g)V2 + d 



n—q 



exp 



2nd 



2 + d 



(i) . 2nd 
Ln - l[ 2+~d> 



(93) 



(94) 



where Eq. ([93]) defines l/n 1 ^ 2 used in appendix C 



Specifically, we employ formula 13.3.8 from 79 with h = 1/2. 
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Finally we give the result for the matrix element (j45| for R = 0. Using the single electron states of the metal 
specified in (142p and measuring length again in units of as, 



I^},(0) — 2 I dxexp[—k s x]sin(kx)sin(k'x) (95) 



= 16k 2 s J {2) (E k ,E k ,)kk' , (97) 
which also defines the function j( 2 \E k , £*')> with £^ = fc 2 and likewise for Ey , needed in appendix C. 

C Calculation of s e and r e 

In this appendix we give mathematical details concerning the calculation of s e and r e . In all equations below we 
use dimensionless variables measuring energies and lengths in units of Rq and ag, respectively. We are furthermore 
interested in the limit L — > oo and A — > oo. Thus, momentum sums become integrals according to 

1 V- f dk , 1 f dQ 



L ^— ' / 2tt A ^ / (2tt) 2 

For the purpose of doing some of the integrals analytically, we found it convenient to rewrite the 6— function for 
energy conservation as follows: 

/oo 
duj5(E Q , q , - E Qn - u)8{E K >k> - E Kk + u) . (99) 
-OO 

The angles can then be integrated out and the global sticking coefficient s e defined in Eq. ([5T]) becomes 



A /q3/2 r oc r oo poo poo poo I r(-0 fri\ r(^) 

i^W dq> dk dk' du dR l J -' (0)J **' 

7T 2 /?e _ JO JO JO J-oo JO 



,2 



2 /?y 2 „ JO JO JO J-oo J + R 2 

x [H-n B (w)]AT(ii,w ) £fc > £; fc /)ii- 1 exp[-/3e*»(ii,w,^)] , (100) 

where, for simplicity, we have neglect the dependence of the electronic matrix elements (|44p and (143)) on the lateral 
momentum transfer R — Q — Q' and introduced two functions: 

N(R,cj,E k ,E k ,) = F_ 1/2 ((3 S {E F - y kk >(R,u) + u) - F_ 1/2 {[3 S {E F - y kk ,(R,u))) , (101) 

= E n +u+ ( E q --E n + R*-u V ( (1Q2) 



2i? 

with i 7, _ 1 / 2 (ic) Fermi integrals for which, as far as the numerics is concerned, we take Unger's approximation [81j . and 

y w (R, U ) = E k + ( Ek ' Ek 2R R2 ' UJ ) 2 (103) 

with E k = k 2 and likewise for E k and E q i . The functions (0) and Iujj, (0) are, respectively, defined in Eqs. ([57)1 
and (|95[) in appendix B. 

Using Eg/, E k , and E^/ instead of q 1 , k, and k' as integration variables, we finally find the result presented in 
Eqs. (H2D and (53]) with 

/•OO 

- 51 \li 1} \ 2 eM-0e(E n +u>)} / dxexp[-/3 ea; 2 ] , (104) 

n J X n {R,Uj) 

p OO 

g(R,uj)= dEdE'J^{E,E')N{R,cu,E,E r ), (105) 
Jo 



and 



*n(i^) = - § - ■ (106) 
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The calculation of the energy resolved sticking coefficient proceeds along the same lines. For perpendicular incidence 
we find the result stated in Eq. (l54|) of the main text with 

gHR, E>) = £ \^\ 2 l + n ff {R / k y R2) 9(R, E' — E n — R 2 ) . (107) 

Now we turn our attention to the calculation of the desorption time. It is quite similar. An intermediate expression, 
after expressing energy conservation in the form (|99p and performing the integrals over angles, is 



8^HZ^fJ Q Vo Jo J-oo Jo & S +B? 
x [l + n B (w)]iV'(i2,w,Sfc,Sfc')i2~ 1 exp[-A#„/(iJ,w,^,)] , (108) 

where we have again neglected the dependence of the electronic matrix elements ([44)) and (|45|) on the transfer of 
lateral momentum and introduced 

Z = Y J cM-PsE n ] , (109) 

$AR,",E q ) =E q + u+( ^ En 2R R2 +^ 2 . (110) 

Using again Ek = k 2 7 Ey = k' 2 , and E q — q 2 as integration variables we finally find the result (|58|) and ([59]) given 
in main text with 

/■oo 

/(£,u;) = V|lW|% X p[-/3 e £ n ] / dyexph^y 2 ] (111) 



and 



y n (R,uj) = — . (112) 



At the end of this appendix let us say a few words about Laplace's approximation [82j which we used to derive 
Eqs. (|6ip - (|63p . If there is only a single bound state with energy E\ the summations over n reduce to a single term. 
For ksTs <C \E\\ and ksT e <C \E\\ it is then possible to do some of the integrals defining s e and r e asymptotically 
within Laplace's approximation. 

First, we consider Laplace's approximation for r e . For a single bound state 

y* oo 

f(R,u) = \4 1] \ 2 expl-PeEx] / dyexp[-(3 s y 2 } . (113) 

Provided ksT s <C \E\\, f(R,ui) is largest for y 1 (R,cu) < 0, that is, for to < E\ — R 2 < 0. In this domain, Laplace's 
approximation to the y— integral gives y7r/ / 9 s /2, where the factor 1/2 anticipates that the R- and cj-integrations are 
later performed also within Laplace's approximation. Changing lo — > — lo, we obtain 

tfcsorb « ^ /" ^ 1 ,7i^ xa lJ^r^Ph^t-g^-^)] > ( 114 ) 

2 Pi' Jo J^i+m 1+UvM 



where we used 1 + w) = — ub(lo). Hence, using ([58 



r 47r 5 / 2 /3l /2 ft 
^ « - > (H5) 



where 



.71 £\|> = / ,/7? / do; . (116) 
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Since ksT s <C we can approximate in Eq. (|116|) the Bose distribution function by exp[— /3 s u>]. Hence, the 

main contribution to the to— integral will come from its lower boundary. Calculating the lo— integral within Laplace's 
approximation and then applying, in a last step, Laplace's approximation also to the remaining R— integral, we find 

JdEiD^^geM-ME^} (117) 

with 

g= Urn -|£ x |)] , (118) 

which, combined with Eq. (|115p . leads to Eq. (|62[) given in the main text. 

Calculating s e within Laplace's approximation is quite similar. However, whereas for r e it is a reasonable ap- 
proximation, because ksT s <C \E\\, for s e it is only meaningful when ksT e is also much smaller than \E\\. Under 
this assumption, which is of course usually not satisfied, we find from Eq. (|53|) . again anticipating that the R- and 
w-integrations are later performed within Laplace's approximation, 

I^^\^\l { i ] \ 2 ^m\Ei\]K{\Ei\) (119) 

Pe 

with 

K(\Et\)= dR te TZTwk L eM-Mg(R,u) , (120) 

Jo JiE^+Ri i + {R/k s y 

to which we again successively apply Laplace's approximation to find 

1/2 

KdEtDn—^g-eM-PelE^} (121) 
2/3 e 3/2 



with 



g= Urn [<?(#, I^il)] . (122) 



Hence, combining (|12ip and (|1 19[) and inserting the result in (|52p gives 



L 4|/W| 2 g 

Se " n0Ve /2 ' (123) 



Using the properties of the function g(R,uj) we now show that g = g. First, we see from the definition ()105jl that 
the R— dependence of g(R, \E\\) comes from the R— dependence of the function ykk'(R, \Ei\) defined in (|103p . Then 
we notice that 

( coforE^E'+|Ei| 
\imy kk ,(R,\E 1 \) = \ , (124) 

[ E for E = E' + |Ei| 

from which follows 

forE^E' + |Ei| 

lim N(R, \Ei\,E,E') = { , (125) 

F- 1/2 (0.(E F - E+ [Et])) - F_ 1/2 (ME F - E)) for E = E' + |Ei| 



R^O 



because F_i/ 2 (x) vanishes for x — > — oo, and thus 

g=limg(R,\E 1 \) (126) 

= / dE 1 J^iE' + \E l lE')[F_ l/2 {(3 s (E F - E 1 )) - F_ l/2 {13 S {E F - E' -\E 1 \))} (127) 
Jo 

poo 

= / dE'jW(E\E' + \E 1 \)[F- 1/2 (/3 a (E F -E'))-F- 1/2 (p a {E F -E'-\E 1 \)j\ , (128) 
Jo 
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where in the last line we used jW(E,E') = JW(E',E). To calculate g we proceed in the same way, noticing however 
that N(R,-\Ex\,E, E') is finite only for E' = E+\E 1 \. Hence, 

g=lim[-g(R,-\E 1 \)} (129) 

xX — >U 

= - / dEjW(E,E + \Ei\)[F- 1/ 2(J3 a (E F -E-\E 1 \))-F- 1/2 (l3 a (E F -E))\ (130) 
Jo 

= dEjW(E,E + \Ei\)[F_ 1/2 (0 a (E F -Er))-F_ 1/2 (J3 a (E F -E-\E 1 \)j\ (131) 
Jo 

= 9 ■ (132) 
Since g — g, Eq. (|123p is identical to Eq. ([6Tj) given in the main text. 
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